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EXECUTIVE  SUMMARY 

We  report  on  the  technical  progress  in  the  first  year  of  an  interdisciplinary  experimental,  numeri¬ 
cal  and  theoretical  program  to  extend  the  state  of  knowledge  and  understanding  of  the  effects  of 
chemical  reactions  in  hypervelocity  flows.  The  program  combines  unique  experimental  and  compu¬ 
tational  facilities  at  Caltech  to  address  the  key  problems  in  aerothermochemistry  that  arise  from 
interactions  between  the  three  strongly  nonlinear  efFects:  Compressibility;  vorticity;  and  chemistry. 
Fig.  1  presents  an  overview  of  the  program. 


Project  /  Investigator 

2 

I 

d 

% 

i 

1 

5 

s 

1 

5 

i 

1 

2 

f 

1.1.  TraiuBlonnl  and  TurbuM  Boundary  Layon 

□ 

□ 

□ 

□ 

□ 

• 

1.2.  NonequBbtfum  and  VorfieHy  Downrireom  of  Bow  Shock, 

□ 

8 

□ 

□ 

■ 

□ 

□ 

□ 

LS.  Mahon  WMi  Chemicol  Reaction 

□ 

□ 

8 

□ 

□ 

r 

• 

□ 

1.1.  NonoquBbrfum  Chembtry  In  Shock-Vortex  Interaction 

■ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

n 

jiwti  "»»w  trmrocnom  n  fiyporriioc/ijr  « tow 

• 

i 

□ 

□ 

0 

□ 

□ 

M.  Superior*:  Shear-Layer  Flow, 

• 

"1 

• 

fii 

■ 

□ 

IV.  1.  Moaturamont  of  Vibrational  Nonequflbrium  Reaction  Rato 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

IV.2.  Rato,  ol  Ooctran  Oman  Roacfion,  In  Hyporaonlc  Flow 

□ 

□ 

L 

■ 

□ 

□ 

NJ.  NonequBbrtum  looward  Shock-Vortex  Aerodynamic, 

■ 

□ 

□ 

□ 

□ 

□ 

8 

V.2.  Until  for  Compaction  and  Temperature  ttooxuremont 

□ 

□ 

□ 

_ 

□ 

□ 

□ 

_ 

■ 

* 

Fig.  1.  Sub-projects,  principal  investigators,  and  main  interactions. 


The  program  has  led  to  the  following  important  new  results  (partial  list): 

•  Extensive  experimental  data  on  the  effect  of  enthalpy  on  transition  in  hypervelocity  flow.  First 
visualization  of  high-enthalpy  turbulent  boundary  layer  structure. 

•  Extensive  new  data  on  hypervelocity  flow  over  spheres,  confirming  numerical  and  new  theoret¬ 
ical  results  as  well  as  T5  enthalpy. 

•  Progress  on  methods  of  parallel  computation  of  shock-vortex  interaction.  New  computations 
of  3-D  leeward  reacting  flow. 

•  Development  and  tests  of  a  greatly  improved  low-cost  thermocouple  heat  transfer  gauge. 

•  Upgrade  of  the  supersonic  shear-layer  facility  to  all-hyperbolic  operation.  Experimental  explo¬ 
ration  of  this  regime  with  a  new  Rayleigh  scattering  method. 

•  Test  of  the  validity  of  vibration-dissociation  coupling  models  from  T5  data  and  computation. 

•  Successful  application  of  new  parallelized  algorithm  of  the  extended  Schwinger  multi-channel 
method  to  computation  of  electron-molecule  collision  cross-sections. 

•  Prediction  and  realization  of  a  new  non-intrusive  diagnostic  method,  laser-induced  thermal 
acoustics  (LITA),  for  accurate  measurement  of  sound  speed  and  bulk  viscosity. 

The  research  has  benefited  extensively  from  close  interaction  within  the  program,  including  regu¬ 
lar  research  conferences  and  seminars,  and  through  extensive  communication  with  outside  groups 
through  presentations  at  conferences  and  seminars,  as  well  as  through  consultancies.  In  most  of  the 
sub-projects  the  progress  made  at  least  measured  up  to  our  initial  plans,  and  in  some  cases  went 
well  beyond,  through  favorable  developments  during  the  year.  The  research  has  led  to  3  Ph.D. 
degrees  and  produced  14  publications.  The  sub-projects  are  actively  producing  results  at  present. 
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INTRODUCTION 


Background 

At  the  outset  of  this  report  it  is  pertinent  to  quote  three  of  the  opening  paragraphs  of  the 
proposal  that  led  to  the  current  URI  Grant: 

“Transport  to  or  from  space,  aero-assisted  orbital  transfer  operations  and  other  hypervelocity 
flight  maneuvers  necessarily  occur  at  such  high  (near  orbital)  speeds,  that  the  bow-shock  heating  of 
the  components  of  air  causes  them  to  be  dissociated,  to  react  with  each  other,  and,  at  particularly 
high  speeds,  to  be  ionized.  Efficient  air-breathing  propulsion  of  such  high-speed  vehicles  also 
requires  supersonic  combustion  of  the  fuel.  To  a  large  degree,  these  applications  have  motivated 
the  interest  in  the  field  of  study  that  has  been  named  aerothermochemistry  by  von  Karman. 

The  character  of  the  dynamics  of  reacting  gases  is  changed  relative  to  perfect-gas  dynamics  in 
two  important  ways:  First,  the  equilibrium  behavior  is  much  more  complicated,  and  second,  charac¬ 
teristic  times  are  introduced  by  the  fact  that  the  reactions  proceed  at  finite  rates.  These  two  types 
of  highly  nonlinear  gas  imperfections,  interacting  with  the  classical  nonlinear  manifestations  of  gas- 
dynamics,  namely  shock  waves  and  turbulence,  produce  many  new  flow  structures  and  phenomena. 
Some  of  these  are  known  and  qualitatively  understood,  a  few  are  quantitatively  predictable  by 
computation,  all  are  of  scientific,  and  most  of  practical  importance  to  the  operations  listed  above. 
The  experimental  data  in  this  field  Me  very  limited,  because  of  the  difficulties  associated  with 
laboratory  simulation.  This  constitutes  the  most  important  need  in  aerothermochemistry. 

The  proposed  research  work  may  be  grouped  under  the  following  major  headings: 
supersonic  turbulent  mixing  and  combustion 
effects  of  chemical  reactions  in  external  hypervelocity  flows 
reaction  rate  models 

The  philosophy  of  this  proposal  is  to  concentrate  interdisciplinary  effort  on  particular  problems  in 
the  first  two  categories.  The  main  thrust  of  the  proposed  research  is  to  use  the  two  new  facilities  (T5 
hypervelocity  shock  tunnel  and  S3LSupersonic  Shear-Layer  Facility)  for  an  experimental  program 
to  study  and  understand  the  effects  introduced  by  chemical  reactions.  A  parallel  theoretical  and 
computational  effort  will  help  design  the  experiments  and  test  models  and  mechanisms  of  the  gas 
behavior.  The  third  category  pervades  the  whole  program  and  represents  research  in  theoretical 
and  computational  chemistry  and  fluid  dynamics,  as  well  as  comparison  with  experiment.” 

The  program  was  then  subdivided  into  the  following  four  parts 

I.  Reaction-rate  controlled  shear  flows 

II.  Shock-vortidty  interactions 

III  Supersonic  mixing  and  combustion 

IV  Chemistry  in  non-uniform  flow, 

all  of  which  involve  complicated  and  non-linear  interaction  between  vorticity,  chemistry  and  com¬ 
pressibility.  By  focusing  the  available  facilities  and  expertise  on  the  main  problems  in  the  field, 
sub-projects  were  determined  as  shown  in  Fig.  refxint  together  with  the  Pi’s  and  the  main  inter¬ 
actions  between  them. 
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A  brief  outline  of  the  objectives  was  then  formnlated  in  the  following  statement  of  work.  The 

numbering  of  these  tub-projects  corresponds  to  that  of  the  sections  in  this  report  under  which  the 

progress  is  described  in  more  detail. 

Statement  of  Work 

1.1  Obtain  qualitative  information  on  the  structure  of  transitional  and  turbulent  boundary  layers 
at  high  enthalpy.  Determine  the  response  of  a  laminar  boundary  layer  to  disturbances  at  the 
most  strongly  amplified  frequency  for  the  second  mode  instability.  Perform  T5  experiments  on 
a  sharp  slender  cone,  using  differential  interferometry,  laser-induced  fluorescence  and  surface 
instruments. 

1.2  Design  experiments  for  testing  the  limits  of  binary  scaling  by  making  a  numerical  study  of 
blunt  body  flows  in  nitrogen  and  air.  Carry  out  T5  experiments  on  spheres  using  surface 
instrumentation,  differential  interferometry  and  degenerate  four-wave  mixing  (DFWM).  Extend 
the  experiments  to  study  the  high-vorticity  layer  downstream  of  a  bow  shock  associated  with 
the  large  density  rise  caused  by  dissociation. 

1.3  Determine  the  flux  of  reactive  and  inert  species,  and  of  heat  and  momentum  to  the  surface 
of  ablating  bodies  in  hypervelocity  flow.  Carry  out  T5  experiments  using  sampling  mass- 
spectrometric  analysis  and  surface  heat  transfer  and  pressure  measurements.  Formulate  ana¬ 
lytical  models  of  reactive  effects  on  ablation. 

II.  1  Perform  ultra-high  resolution  numerical  experiments  combined  with  nonlinear  analysis  and 
numerics  to  establish  appropriate  initial  conditions  to  elucidate  the  generic  mechanisms  of 
shock- vortex  interactions,  namely,  (1)  interaction  of  a  shock  with  a  vorticity  field  in  the  presence 
of  nonequilibrium  chemistry,  and  (2)  vorticity  generation  in  shock  focusing  with  chemistry. 

II.2  Determine  the  effects  of  chemical  relaxation  on  the  shock-on-shock  problem.  Demonstrate  the 
magnitude  of  nonequilibrium  effects  on  shock  impingement  heating.  Carry  out  experiments 
with  a  2-dimensional  model  in  T5,  supported  by  local  shock-polar  analysis  and  numerical 
calculations  of  an  ideal  dissociating  gas. 

Ill  Carry  out  S3L  experiments  with  chemically  reacting  and  non-reacting  flows,  using  schlieren, 
spectroscopy,  Rayleigh  scattering,  chemical  effects  and  new  diagnostic  methods,  to  gain  under¬ 
standing  of  the  mechanisms  of  mixing  when  a  shock  interacts  with  a  supersonic  turbulent  shear 
layer,  and  in  the  flow  downstream  of  an  expansion  wave.  Develop  new  related  diagnostics  and 
computational  capabilities. 

IV.l  Design  and  carry  out  an  experiment  in  T5  to  provide  data  to  test  chemical  reaction  models  for 
flows  with  vibrational  nonequilibrium.  Use  a  previously  developed  computational  technique 
to  simulate  the  flow  over  proposed  experimental  configurations.  Identify  a  configuration  that 
exhibits  significant  effects  of  vibrational  nonequilibrium  on  the  reaction  rates. 

IV.2  Exploit  computational  methodologies  and  high-performance  and  cost-effective  computing  pro¬ 
vided  by  highly  parallel  and  scalable  supercomputers  to  solve  the  complex  equations  which 
govern  electron-molecule  collision  processes  needed  for  robust  modeling  of  the  chemical  and 
physical  properties  of  hypervelocity  flows. 

FV.3  Demonstrate  the  effect  of  chemistry  on  the  leeward  shock-vortex  separation  region  of  hyper- 
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velocity  flow  about  slender  bodies  using  high-resolution  CFD  on  the  Intel  Touchstone- Delta 
computer,  combined  with  dimensional  analysis  and  physical  scaling. 

V.2  Develop  DFWM  for  measuring  species  concentration  and  temperature  in  hypervelocity  flows, 
concentrate  on  NO-concentrations  because  of  their  importance  in  air  flows.  Buy  the  remaining 
equipment  to  complete  the  set  of  components,  assemble  the  instrumentation,  debug  and  apply 
it  to  T5  flows. 

Current  Status  of  Progress 

In  the  following  paragraphs,  which  mirror  the  statement  of  work,  the  progress  in  the  sub- 

projects  is  summarized. 

1.1  Extensive  experiments  were  conducted  over  the  enthalpy  range  3.5MJ/kg  <  /»o  <  22MJ/kg  in 
air  and  nitrogen,  and  3.5  MJ/kg  <  ho  <  8.6MJ/kg  in  carbon  dioxide,  with  laminar,  transitional 
and  turbulent  flow  on  a  slender  cone.  Flow  visualization  with  a  novel  method  and  heat  flux 
measurements  show  that  the  mechanism  of  transition  under  the  conditions  examined  is  via 
the  Tollmien-Schlichting  instability.  High-enthalpy  real-gas  effects  have  a  strongly  stabilizing 
effect,  which  appears  to  scale  with  total  enthalpy  normalized  by  dissociation  energy. 

1.2  A  series  of  experiments  on  flow  over  spheres  has  been  completed  that  covers  the  full  range  of 
enthalpy  of  T5  in  air,  nitrogen  and  COo  using  interferometry  and  heat  flux  measurements.  Nu¬ 
merically  computed  interferograms  agred  well  with  experiment  except  in  small  regions  near  the 
shock  (see  also  IV. 1).  Theoretical  analysis  of  the  problem  led  to  a  complete  understanding  of 
the  correlation  of  the  dimensionless  stand-off  distance  with  dimensionless  reaction-rate  param¬ 
eter,  and  its  extension  to  reacting  gas  mixtures.  A  previous  theoretical  analysis  of  dissociating 
flow  along  a  streamline  downstream  of  a  curved  shock  has  been  extended  to  include  the  effects 
of  recombination.  The  limits  of  binary  scaling  were  delineated.  Initial  exploratory  experiments 
on  the  detection  of  the  high  vorticity  layer  and  its  stability  were  performed. 

1.3  This  subproject  is  not  scheduled  to  start  until  the  end  of  the  second  year,  when  it  is  planned 
that  it  will  follow  the  completion  of  the  work  under  II.2.  preliminary  plans  for  the  mass- 
spectrometer  design  have  been  made  but  are  not  reported  in  detail. 

II.  1  Progress  has  been  made  in  computation  of  steady  compressible  vortex  structures  which  will 
then  be  used  as  special  initial  conditions  for  controlled  numerical  experiments  on  shock  vor¬ 
ticity  interaction.  Progress  has  also  been  made  in  developing  the  computational  capability 
for  simulation  of  vorticity  generation  through  shock  focusing.  A  large  scale  two  dimensional 
compressible  flow  solver  has  been  constructed  and  a  high  resolution  study  has  been  started 
on  two  massively  parallel  machines  including  the  Intel  Delta  and  the  IBM  SP/1.  Whitham 
theory  was  applied  to  the  computation  of  an  approximation  for  the  vorticity  generated  by  a 
shock  as  it  interacts  with  a  density  gradient.  Computations  are  now  in  progress  to  test  this 
approximation. 

II.2  The  primary  effort  during  this  reporting  period  was  to  make  major  improvements  in  methods 
for  heat  transfer  measurement  and  optical  interferometry  in  T5.  In  addressing  the  planned 
improvement  of  the  existing  Wollaston  prism  differential  interferometer,  a  more  sweeping  solu¬ 
tion  was  sought;  a  completely  new  holographic  system  for  optical  diagnostics  in  T5  has  been 
designed  and  built.  This  instrument  will  be  reported  at  another  time.  In  this  report  we  doc- 
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ument  the  development  of  a  new  heat  transfer  gauge  which  gives  performance  superior  to  any 
existing  device. 

Ill  In  the  compressible  shear-layer  turbulence  and  supersonic  mixing  part  of  this  effort  we  have 
upgraded  the  supersonic  shear  layer  (S3L)  facility  to  accommodate  bi-supersonic  flow,  i.e., 
both  freestreams  supersonic,  resulting  in  an  all-hyperbolic  flow  configuration.  A  first  set  of 
experiments  were  performed  in  this  flow  regime.  A  second  generation,  higher  signal-to-noise 
ratio,  planar,  laser  Rayleigh-scattering  experiments  were  performed  in  supersonic  shear  layers 
with  varying  degrees  of  compressibility,  which  indicate  the  existence  of  a  high-compressibility 
flow  regime  with  different  entrainment  and  mixing  mechanisms.  Computationally,  we  have 
performed  simulations  of  bi-supersonic  and  transonic  shear-layer  flows  that  reproduce  many  of 
the  features  observed  in  our  preliminary  experiments  in  this  flow  regime. 

IV.l  All  currently-used  vibration-dissociation  coupling  models  have  been  implemented  in  a  two- 
dimensional  thermo-chemical  nonequilibrium  computational  fluid  dynamics  method.  The  flow 
of  reacting  nitrogen  over  a  sphere  has  been  simulated  using  these  models  and  compared  to 
available  experimental  data  from  T5.  Under  the  conditions  of  the  flow  tested,  it  was  found 
that  while  the  various  models  gave  different  vibrational  temperature  distributions,  the  shock 
stand-off  distances  and  density  contours  are  very  similar.  A  new  experiment  is  being  designed 
for  T5  that  will  be  more  sensitive  to  the  differences  in  the  reaction  rate  models.  Simulations 
have  been  made  of  T  5  flows  over  spheres  at  lower  densities,  and  differences  among  the  models 
are  discernible  in  both  the  shock  stand-o,f  distances  and  the  density  contours. 

IV.2  The  computationally  most  intense  step  of  the  previously  developed  program  for  computation 
of  low-energy  electron- molecule  collisions  with  the  extender  Schwinger  multi-channel  method 
has  been  parallelized.  Sustained  performance  has  been  achieved  at  about  4  GFLOP  on  the  full 
Intel  Delta,  approximately  100  times  the  performance  of  the  original  sequential  program  on 
the  CRAY  Y-MP.  The  program  has  been  used  in  the  computation  of  a  number  of  important 
collision  processes  and  cross-sections. 

V.2  In  this  sub-project,  a  new  diagnostic  method,  laser-induced  thermal  acoustics  (LITA),  was 
discovered  which  is  more  powerful  than  the  originally  planned  DFWM.  The  theory  of  LITA 
was  developed,  and  this  predicted  that  the  technique  would  permit  non-intrusive  measurement 
of  the  speed  of  sound  and  the  bulk  viscosity  of  a  gas.  Extensive  equipment  was  bought  and 
the  method  was  developed  and  tested.  The  technique  proved  to  be  entirely  successful  and 
very  accurate,  with  many  advantages  over  other  laser  diagnostic  techniques,  and  extensive 
further  possibilities.  The  theory,  developed  six  months  before  the  tests,  precisely  predicted  the 
recorded  signals  in  the  tests. 

A  total  of  4  higher  degrees  were  earned.  3  papers  were  published  or  accepted  for  publication  in 
the  archival  literature,  and  5  have  been  submitted  or  are  in  preparation.  10  papers  were  published 
as  conference  papers  or  in  proceedings.  The  work  has  resulted  in  a  large  number  of  presentations 
at  various  seminars,  meetings  and  conferences.  Details  of  these  are  given  in  the  subsections. 

Interactions 

The  interaction  between  the  principal  investigators  and  their  research  groups  within  the  project 
followed  the  diagram  in  Fig.  1.  very  closely.  Much  of  this  interaction  is  forced  by  joint  use  of 
facilities  or  instrumentation,  sharing  of  technical  personnel,  or  through  discussions  that  arise  at 


almost  weekly  interactions  in  the  Fluid  Mechanics  Research  Conference  at  GALCIT.  Interaction 
between  the  separate  Aerot hemodynamics  UltTs  benefited  from  the  meeting  held  at  Caltech  in 
July  1993.  Advantage  was  taken  by  several  Pi’s  of  attendance  at  conferences  (AIAA,  19thISSW) 
and  meetings  at  AEDC  to  interact  with  other  groups  nationally  and  internationally. 

Layout 

The  following  chapters  present  the  main  results  of  the  work  in  more  detail,  in  the  order  shown 
above.  Each  subproject  then  lists  the  personnel,  publications,  degrees  earned,  and  interactions 
separately.  A  few  of  the  subprojects  had  received  preliminary  funding  before  the  URI  funding 
began,  so  that  they  are  at  a  more  mature  stage  than  others. 


CHAPTER  I 

REACTION-RATE  CONTROLLED  SHEAR  FLOW 


1.1  TRANSITIONAL  AND  TURBULENT  BOUNDARY  LAYERS 

Objectives  and  Status  of  Research 

The  objective  of  this  experimental  investigation  is  to  obtain  qualitative  information  on  the 
structure  of  transitional  and  turbulent  boundary  layers  at  high  enthalpy,  and  to  determine  the 
response  of  the  laminar  boundary  layer  to  disturbances  at  the  most  strongly  amplified  frequency 
for  the  second  mode  instability.  The  experiments  use  surface  heat  transfer  measurements  on  a 
slender,  sharp  cone  and  optical  flow  visualization  to  detect  effects  of  chemical  reactions  on  the  heat 
transfer,  transition  location,  and  structure  of  the  flow  field.  A  gas-dynamical  oscillator  is  to  be 
used  to  generate  the  controlled  disturbances. 

The  proposal  for  this  sub-project  of  the  grant  was  written  on  the  basis  of  preliminary  experi¬ 
ments  performed  with  the  slender  cone  in  T5  that  were  reported  in  Ref.  1.  During  the  reporting 
period  an  extensive  investigation  covered  the  enthalpy  range  3.5  MJ/kg  <  hQ  <  22MJ/kg  in  air  and 
nitrogen,  and  3.5  MJ/kg  <  ho  <  8.6M J/kg  in  carbon  dioxide.  The  results  include  extensive  data 
on  surface  heat  transfer  on  the  slender  cone  in  laminar,  transitional  and  turbulent  boundary  layers, 
from  the  perfect-gas  to  the  high-enthalpy  regime,  in  the  latter  with  and  without  recombination  at 
the  wall.  A  novel  method  is  employed  to  visualize  the  boundary  layer  flow  structure,  yielding  the 
first  pictures  of  turbulent  boundary  layer  structure  in  high-enthalpy  flow.  The  results  demonstrate 
that  the  similarity  theories  predict  the  heat  transfer  rate  correctly  in  the  laminar  boundary  layer. 
The  mechanism  of  transition  under  the  conditions  examined  is  via  the  Tollmien-Schlichting  insta¬ 
bility.  High-enthalpy  real-gas  effects  have  a  strongly  stabilizing  effect,  which  appears  to  scale  with 
total  enthalpy  normalized  by  dissociation  energy. 

New  experiments  are  now  under  way  to  extend  and  consolidate  these  findings,  and  to  con¬ 
tinue  with  the  much  more  ambitious  stability  investigation,  as  well  as  with  a  more  sophisticated 
visualization  technique. 

1.1.1  The  experimental  setup 

The  tests  were  conducted  in  the  hypervelocity  shock  tunnel  T5,  which  has  been  adequately 
described  in  the  proposal  for  this  URI  grant.  The  slender  cone  used  in  the  experiments  is  1  m 
long  and  is  divided  into  three  sections  as  shown  in  Fig.  1.1  The  tip  section  (76  mm  long)  is  a  solid 
piece  of  Waspalloy,  a  stainless  steel  that  has  good  toughness  at  high  temperature.  Even  so,  the  tip 
radius  was  found  to  increase  progressively  by  melting  from  shot  to  shot,  from  an  initial  0.06  mm  to 
approximately  0.15  mm,  at  which  value  it  stabilizes  after  typically  20  shots.  This  increase  of  nose 
radius  was  found  to  have  no  detectable  effect  on  the  heat  transfer  distribution. 


I 

Fig.  1.1.  Photograph  of  the  model  partly  disassembled  to  show  the  tip  (a),  the  second  section  (b) 
and  the  main  body  (c).  The  sting  mount  is  indicated  as  (d)  and  the  copper  pipe  protecting 
the  fifteen  thermocouple  wires  as  (e). 

The  second  section  (127  mm  long)  is  a  solid  piece  of  aluminum.  The  third  piece  (790  mm  long, 
made  of  aluminum)  is  hollow  and  has  a  wall  thickness  of  12  mm.  It  consists  of  two  longitudinal 
halves  held  together  by  a  dove-tail  joint,  so  that  they  slide  into  each  other  without  any  need  for 
screws,  leaving  the  outside  perfectly  smooth  except  for  an  almost  invisible  longitudinal  seam.  The 
assembled  cone  is  screwed  onto  a  sting-mounted  base  plate  made  of  steel. 

This  is  mounted  rigidly  in  the  test  section  of  T5  after  careful  alignment  of  the  cone  and  nozzle 
axes.  The  model  is  positioned  in  the  test  section  such  that  it  protrudes  330  mm  into  the  nozzle. 
This  allows  all  the  gauges  on  the  cone  to  be  completely  inside  the  uniform  test  rhombus  of  the 
nozzle  flow  above  Moo  =  5.6. 


1. 1.1.1  Heat  transfer  gauges 

Fifteen  thermocouples  of  type  E  (chromel-constantan),  manufactured  by  Medtherm  Corp., 
were  mounted  flush  with  the  surface  of  the  cone  to  measure  the  heat  transfer  distribution.  The 
gauges  are  placed  at  approximately  logarithmic  spacing  in  axial  distance  from  the  tip  from  241  mm 
to  889  mm.  Their  azimuthal  position  varies  in  order  to  avoid  mutual  influences,  though  their 
flushness  is  very  good.  The  temperature  signal  from  the  gauges  is  converted  to  heat  transfer  by 
the  method  of  Schultz  and  Jones  (1973).  For  more  detail  of  the  theory  of  heat  flux  measurement 
see  section  II.2,  where  the  development  of  improved  heat  flux  gauges  is  presented. 
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1. 1.1.8  Flow  visualization 

The  main  interest  in  applying  flow  visualization  to  this  experiment  is  to  study  the  structure  of 
the  transitional  and  turbulent  boundary  layer  in  the  high-enthalpy  regime.  The  density  and  optical 
path  length  are  too  small  for  regular  schlieren,  shadowgraphy  or  interferometry  to  be  sufficiently 
sensitive  to  detect  the  differences  in  density  that  occur  in  the  boundary  layer  because  of  structure. 
We  therefore  resorted  to  the  idea  of  resonant  enhancement  of  the  refractive  index  of  the  medium 
used  by  Blenstrup  et  al.  (1979).  This  technique  relies  on  the  fact  that  the  refractive  index  of  a 
gas  reaches  large  extrema  near  a  spectral  line  of  the  medium.  In  the  visible  range,  the  technique 
may  be  implemented  by  seeding  the  flow  with  sodium  and  using  a  light  source  with  a  wavelength 
just  off  the  center  of  one  of  the  D-lines.  Even  if  the  light  source  wavelength  is  not  perfectly  tuned, 
such  that  it  is  in  the  center  of  the  D-line,  where  the  refractive  index  contribution  of  the  line  is 
zero,  the  absorption  of  the  light  by  the  resonance  is  very  high.  The  technique  will  then  still  work 
satisfactorily  for  qualitative  flow  visualization  with  good  repeatability. 

In  our  experiment,  a  preferential  enhancement  of  the  sensitivity  of  the  boundary-layer  fluid 
was  desirable,  so  we  seeded  the  flow  only  within  the  boundary  layer,  by  depositing  a  very  small 
drop  of  a  saturated  saline  solution  on  the  portion  of  the  cone  surface  of  interest.  This  was  then 
left  to  dry,  so  that  a  thin  film  of  salt  was  left  to  be  ablated  during  the  test.  When  dry,  the  film 
actually  consists  of  cubic  crystals  100  pm  high  spaced  2  or  3  mm  apart. 

The  light  source  used  was  a  tunable  dye-laser  built  for  this  purpose.  It  is  pumped  by  a 
Neodymium:YAG  laser.  Only  one  picture  can  be  taken  during  each  run.  The  laser  is  triggered  by 
an  adjustable  delay  circuit  which  in  turn  is  triggered  by  the  arrival  of  the  shock  at  the  end  of  the 
shock  tube. 

The  resonant  enhancement  of  the  refractive  index  turned  out  to  be  particularly  successful  at 
high  enthalpies  (where  it  was  most  needed  because  of  the  lower  density),  partially  as  a  result  of 
the  relatively  higher  temperature  in  the  boundary  layer,  which  causes  the  salt  to  be  ablated  more 
effectively.  By  control  experiments  it  is  found  that  the  disturbance  caused  by  the  salt  deposition 
is  insufficient  to  cause  the  transition  location  to  be  moved. 

Fig.  1.2  shows  a  schematic  diagram  of  the  model,  the  optical  window,  and  the  limits  of  a  typical 
picture  shown  in  this  report.  This  view  is  always  oriented  parallel  to  the  surface  of  the  cone  and 
always  shows  the  flow  going  from  left  to  right. 


Fig.  1.2.  Schematic  diagram  explaining  the  photographs  to  be  shown  below.  The  flow  is  always  seen 
going  from  left  to  right  and  the  shock  is  out  of  the  field  of  view. 


1.1.2  Heat  transfer  in  the  laminar  regime 
1. 1.2.1  Predictions  of  similarity  theory 


The  Stanton  number  is  defined  as 


St  = 


_ 9w _ 

peuc(hmw  -  hwy 


(i) 


where  qw  is  the  heat  flux  to  the  wall,  pe  and  ue  are  density  and  velocity  at  the  edge  of  the  boundary 
layer,  haw  and  /»«,  are  adiabatic  and  actual  specific  enthalpy  at  the  wall.  To  show  an  example  of 
behavior  differences  between  a  frozen  flow  and  a  chemical  equilibrium  flow,  we  focus  on  the  flow 
condition  corresponding  to  a  specific  enthalpy,  ho,  of  22  MJ/kg  and  stagnation  pressure,  po,  of 
55  MPa.  At  these  conditions,  similarity  theory  gives 


Stair  — 


0.533 


and  St, 


0.536 


Pr2/iRel/2  “lu  J‘"“ro*e"  “  />r2/3*cl/2- 
If  we  include  the  correction  factor  for  real-gas  effects  (chemical  equilibrium), 

0.739  _  0.669 


Stair  = 


PrWRe1/2 


and  St 


nitrogen  — 


pr2/3Rel/2- 


(2) 


(3) 


Re  is  the  Reynolds  number  evaluated  at  the  edge  of  the  boundary  layer  and  Pr  is  the  Prandtl 
number.  We  can  also  compare  these  predictions  with  that  of  low  specific  enthalpy  (same  ior  air 
and  nitrogen)  corresponding  to  ho  =  5.8  MJ/kg  and  to  po  =  15  MPa 


St  = 


0.487 

Pr2/3R€U2- 


Thus,  the  frozen  cases  are  almost  identical. 
1. 1.2.2  Experimented  results 


(4) 


An  example  of  the  normalized  experimental  results  is  shown  in  Fig.  1.3  for  air.  The  data 
shown  as  squares  (po  =  24  MPa  and  ho  «  12  MJ/kg)  indicate  a  frozen  flow  behavior  as  observed 
by  East  et  al.  (1980)  on  a  flat  plate.  The  rms  error  is  of  order  10%.  At  high  specific  enthalpy 
(ho  =  22  MJ/kg),  and  higher  pressure,  the  data  are  more  scattered  (the  rms  error  is  of  order  16%) 
but  fall  clearly  above  the  prediction  of  a  frozen  flow,  indicating  a  higher  heat  release  to  the  wall 
because  of  recombination.  The  presence  of  atomic  species  is  due  to  the  combined  effect  of  free- 
stream  dissociation  and  frictional  heating  in  the  boundary  layer.  The  data  for  air  reach  a  higher 
level  than  corresponding  data  for  nitrogen  because,  at  this  enthalpy,  pure  nitrogen  does  not  get 
dissociated  as  much  as  the  oxygen  in  air,  so  that  there  is  less  recombination  heat  release. 


1.1.3  The  transitional  regime 
1. 1.8.1  Transition  Reynolds 

In  the  example  presented  in  Fig.  1.4  and  Fig.  1.5,  the  normalized  experimental  data  show  that 
the  beginning  of  the  transition  Re%r  corresponds  to  a  Reynolds  number  of  2.5-3  million  whereas  its 
end  corresponds  to  Re  =  4  —  5  million.  For  each  run,  Reu  is  obtained  as  the  intersection  of  two 
lines:  the  first  line  is  a  linear  fit  (of  slope  -0.5)  to  the  first  four  data  points  and  the  second  line  is 
a  linear  least-squares  fit  to  the  data  points  between  this  ‘laminar  line*  and  the  maximum  point. 
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Fig.  1.3.  Normalized  heat  transfer  results  from  two  different  shots  in  T5.  Squares  (run  545):  ho  = 
11.3  MJ/kg  and  po  =  24  MPa;  circles  (run  157):  ho  =  18.5  MJ/kg  and  po  =  60  MPa.  The 

theoretical  predictions  for  these  air  runs  are  also  plotted. - - - :  linear  curve  fit 

of  experimental  data  obtained  from  DiCristina  (1970)  and  from  Owen  et  al.  (1975).  The 
conditions  of  their  experiments  were  ho  <  1  MJ/kg  and  M =  10  and  7,  respectively. 


nitrogen 


Fig.  1.4.  Normalized  heat  transfer  data  from  three  different  shots  in  T5.  Squares  (run  540):  ho  = 
7.7  MJ/kg  and  po  =  57  MPa;  circles  (run  548):  ho  =  9.7  MJ/kg  and  po  =  55  MPa; 
diamonds  (run  534  ):  ho  =  10.6  MJ/kg  and  po  =  55  MPa.  The  theoretical  predictions  for 
run  548  are  also  plotted. 

In  a  given  facility  and  at  given  ambient  noise  conditions,  the  transition  Reynolds  number  may 
be  expected  to  be  related  to  the  following  dimensionless  variables. 

Retr  =  Retr(M,  7,  Pr ,  /?),  (5) 

where  f2  is  a  dimensionless  reaction  rate  parameter,  assuming  there  is  only  one  D  and  only  one  f2. 
For  example,  Demetriades  (1977)  obtains,  for  Tw/Te  =  4.2  and  8.4,  Retr  =  3.3  x  106  and  5.0  x  10®, 
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air 


Fig.  1.5.  Normalized  heat  transfer  measurements  of  three  different  shots  in  T5.  Squares  (run  541): 
ho  =  7.3  MJ/kg  and  po  =  55  MPa;  circles  (run  549):  ho  =  9.0  MJ/kg  and  po  =  52  MPa; 
diamonds  (run  556,  refer  to  run  537  in  Table  2):  ho  =  10.4  MJ/kg  and  po  =  55  MPa.  The 
theoretical  predictions  for  run  549  are  also  plotted. 

respectively.  These  data  and  those  of  DiCristina  will  be  used  for  comparison  purposes1 . 

In  laminar  boundary  layers  of  perfect  gases,  it  is  well  known  that  the  separate  dependence  of 
dimensionless  quantities  on  Af,  7,  Pr  and  Tw /Taw  disappears  if  the  physical  parameters  are  eval¬ 
uated  at  conditions  corresponding  to  the  reference  temperature  T*  introduced  by  Eckert  (1955)2. 
This  concept  worked  very  well  for  the  laminar  regime  on  which  a  dimensionless  correction  factor, 
to  account  for  compressibility,  was  based  (J2ejr  corresponds  to  the  end  of  the  laminar  regime).  T* 
is  given  by 

^  =  0.5  +  0.039M2  +  0.5~.  (6) 

le 


In  an  attempt  to  separate  out  the  real-gas  effects,  i.e.,  the  influence  of  ho/D  and  of  i?,  we 
renormalized  the  results  for  the  transition  Reynolds  number  by  evaluating  it  at  T *  to  give  Re\r. 


Re* 


P*uex 


(7) 


The  results  are  shown  in  Fig.  1.6  together  with  the  data  of  DiChristina  and  Demetriades.  As  may 
be  seen,  our  data  collapse  on  two  lines:  one  for  air  and  one  for  nitrogen. 


Clearly,  the  influence  of  ho  on  Refr  is  very  significant.  Our  data  extends  the  enthalpy  range 
of  tunnel  data  from  a  previous  3  MJ/kg  to  15  MJ/kg.  Also,  the  effect  of  ho  on  the  transition 
Reynolds  number  is  larger  if  the  dissociation  energy  is  smaller.  This  was  further  supported  by  a 
few  experiments  performed  in  carbon  dioxide  (with  a  much  smaller  first  dissociation  energy  and 
large  vibrational  energy),  which  exhibited  a  much  stronger  enthalpy  dependence,  and  with  helium, 
for  which  no  enthalpy  dependence  was  observed. 


1  Not  enough  information  is  provided  by  Martellucd  and  Laganelli  (1974)  to  deduce  Tw/Te-  Demetriades’  (1977) 
data  are  all  obtained  in  the  same  wind  tunnel  and  with  the  same  model. 

2  These  conditions  should  correspond  to  a  reference  enthalpy  h*  rather  than  to  T* .  For  simplicity,  we  assume  that 
the  physical  parameters  p  and  p  evaluated  at  T*  are  good  approximations. 
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Fig.  1.6.  Transition  Reynolds  number  evaluated  at  the  reference  temperature.  Our  data  are  shown 
as  square  symbols  (open  for  air  and  solid  for  nitrogen).  The  error  bar  shown  is  an  ap¬ 
proximation  of  the  accuracy  in  determining  Retr.  The  data  shown  as  open  circles  come 
from  DiCristina  and  Demetriades.  The  two  lines  are  least-squares-fit  approximations  of 
the  open  or  the  solid  symbols  respectively. 

1.1.4  Flow  visualization 

In  this  section  we  show  examples  of  flow  visualization  in  the  laminar,  transitional  and  tur¬ 
bulent  regime.  Fig  1.7  shows  a  digitally  enhanced  small  section  of  the  laminar  boundary  layer, 
visualized  with  the  resonantly  enhanced  sodium  line  technique  at  a  Reynolds  number  of  approxi¬ 
mately  1.2  million.  The  flow  is  clearly  laminar. 


Fig.  1.7.  Flow  of  nitrogen  at  ho  =  12  MJ/kg  and  po  =  24  MPa  (run  544).  This  resonantly  enhanced 
shadowgraph  shows  the  boundary  layer  on  the  top  of  the  5  deg.  half-angle  cone  at  about 
570  mm  from  the  cone  tip  (corresponding  to  Re  =  1.2  million).  The  horizontal  width 
shown  corresponds  to  approximately  12  mm.  In  the  lower  part  of  the  boundary  layer  the 
source  light  is  completely  absorbed. 


In  the  transitional  regime,  it  is  necessary  to  show  a  larger  extent  of  the  boundary  layer.  This  is 
done  in  Fig  1.8.  As  may  be  seen,  the  photograph  exhibits  a  waviness  at  the  left,  which  develops  into 
an  almost  turbulent  boundary  layer  at  the  far  right  of  the  upper  picture.  The  wavelength  (1.5  cm) 


Fig.  1.8.  Flow  of  nitrogen  at.  Iiq  =  11  MJ/kg  and  pa  —  60  MPa  (run  360)  This  resonantly  enhanced 
interferogram  shows  the  boundary  layer  on  the  top  of  the  cone,  starting  on  the  left  from 
x  -  510  mm  (corresponding  to  Re  -  2.75  million)  to  x  =  680  mm  (Re  =  3.77  million). 
The  pictures  below  magnify  portions  of  the  one  above  three  times  (3x).  A  wave  in  the 
laminar  part  of  the  flow  appears  distinctly  and  its  wavelength  is  approximately  1.5  cm. 

corresponds  approximately  to  the  Tollmicn-Schlichting  frequency  at  half  the  free-stream  velocity. 

Finally,  Fig  1.9  shows  a  turbulent  boundary  layer  that  exhibits  clear  oblique  structures  that 
are  suggestive  of  the  kind  seen  in  subsonic  turbulent  boundary  layers.  The  thick  viscous  sublayer 
that  characterizes  cold  hypersonic  flows  is  not  in  evidence  at  these  conditions  at  all.  This  may  be 
because  the  wall-to- reservoir  temperature  ratio  is  always  very  low  in  these  experiments,  or  because 
of  real-gas  effects.  However,  one  also  has  to  be  very  careful  in  the  interpretation  of  the  pictures, 
as  they  are  line-of-sight  integrating  in  an  axisymmctric  flow.  Thus,  the  light  beam  that  passes 
tangentially  to  the  surface  in  the  vertical  ccnterplane  of  the  cone,  also  traverses  portions  of  the 
boundary  layer  at  larger  distances  *Voin  the  wall  before  and  after  traversing  the  ccnterplane. 
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1.2  NONEQUILIBRIUM  AND  VORTICITY  DOWNSTREAM  OF  BOW  SHOCKS 


Objectives  and  Status  of  the  Research 

The  objectives  of  this  sub-project  of  the  program  are  to  design  and  execute  a  series  of  experi¬ 
ments  as  well  as  to  perform  numerical  computations  and  theoretical  analysis,  in  order  to  test  the 
nature  and  limits  of  binary  scaling,  and  to  understand  and  extend  the  reaction-rate  correlation  for 
the  shock-wave  stand-oif  distance,  in  reacting  blunt-body  flows.  In  addition,  the  experiments  are  to 
be  extended  to  investigate  the  high-vorticity  layer  downstream  of  the  bow  shock  that  is  associated 
with  the  large  density  rise  caused  by  dissociation.  The  objectives  have  been  extended  relative  to 
those  in  the  original  proposal  because  of  favorable  developments  during  the  reporting  period. 

An  extensive  series  of  experiments  on  flow  over  spheres  has  been  completed.  This  involved 
nitrogen,  air  and  carbon  dioxide  flows  covering  the  full  range  of  enthalpy  of  the  hypervelocity 
shock  tunnel  T5,  ranging  from  perfect-gas  to  near-equilibrium  flow,  and  a  large  range  of  pressure. 
Differential  interferometry  and  surface  heat  transfer  constituted  the  diagnostics  used.  The  exper¬ 
iments  were  accompanied  by  detailed  numerical  computations  of  the  flow  fields,  which  were  used 
to  construct  interferograms  for  comparison  with  the  experimental  ones.  Theoretical  analysis  of  the 
problem  led  to  a  complete  understanding  of  the  correlation  of  the  dimensionless  stand-off  distance 
with  dimensionless  reaction-rate  parameter,  and  its  extension  to  reacting  gas  mixtures.  Numeri¬ 
cal  and  experimental  results  support  the  correlation  very  well.  A  previous  theoretical  analysis  of 
dissociating  flow  along  a  streamline  downstream  of  a  curved  shock  has  been  extended  to  include 
the  effects  of  recombination.  The  limits  of  binary  scaling  were  delineated.  Initial  exploratory 
experiments  on  the  detection  of  the  high  vorticity  layer  and  its  stability  were  performed. 

Hypervelocity  Flow  over  Spheres 

The  nature  of  the  nonequilibrium  flow  of  dissociating  gases  over  spheres  was  investigated  ex¬ 
perimentally,  numerically  and  theoretically.  A  series  of  experiments  with  three  different  gases, 
nitrogen,  air  and  carbon  dioxide,  was  performed  in  the  shock  tunnel  T5  at  GALCIT.  Five  spheres 
of  different  diameters  (1,  2,  3,  4  and  6  in.)  equipped  with  thermocouples  for  surface  heat  flux 
measurements  were  used.  The  numerical  method  by  Candler  (1988)  was  used  to  conduct  a  parallel 
study  which  strongly  complemented  the  experimental  and  theoretical  efforts. 

1.2.1  Interferometry  and  heat  flux  measurements  of  blunt  body  flows 

The  experimental  differential  interferograms  were  compared  with  the  images  constructed  from 
computational  flow  fields.  Generally,  good  agreement  of  fringe  pattern  and  shock  shape  was  ob¬ 
served.  Fig  1.10  shows  the  comparsion  of  experimental  and  computational  differential  interferograms 
of  the  flow  of  nitrogen  over  a  4  in.  diameter  sphere  at  a  specific  reservoir  enthalpy  of  16  MJ/kg 
and  a  reservoir  pressure  of  60  MPa.  As  may  be  seen,  the  features  of  the  photo  are  faithfully  repro¬ 
duced  by  the  computation,  both  qualitatively  and  quantitatively  in  most  of  the  flow  field.  Small 
differences  near  the  shock  give  valuable  information  about  the  beginning  of  the  dissociation  process. 

Experimental  heat  flux  measurements  show: 

1.  Good  agreement  among  the  measured  stagnation  point  heat  transfer  rates,  computational 

results  and  Fay  and  Riddell’s  theoretical  predictions.  This  also  provides  confirmation  of  T5 

enthalpy 


Fig.  1.10.  Comparison  of  experimental  and  computed  differential  interferograms  of  nitrogen  flow 
over  a  sphere.  Free-stream  conditions:  5.08  km/s,  0.04  kg/in3,  16.56  MJ/kg.  The  pho¬ 
tograph  on  the  right  is  a  finite-fringe  differential  interferogram  of  nitrogen  flow  over  a 
sphere.  The  picture  on  the  left  is  a  corresponding  computed  interferogram  at  the  same 
conditions  as  the  experiment.  Except  for  a  slight  difference  in  the  vicinity  of  the  shock, 
the  two  pictures  are  virtually  congruent.  To  show  this,  a  line  along  the  center  of  the 
calculated  white  fringe  is  superimposed  on  the  photograph. 

2.  For  nitrogen  and  air,  the  measured  heat  flux  distributions  were  also  in  good  agreement  with 
numerical  results  and  Lees’  theory. 

3.  For  carbon  dioxide,  large  differences  between  experimental  and  theoretical  heat  flux  distribu¬ 
tions  were  observed.  Early  transition,  tripped  by  surface  roughness  is  a  possible  cause. 

1.2.2  Binary  scaling  law 

An  analytical  solution  was  obtained  for  inviscid  hypervelocity  dissociating  flow  over  spheres. 
The  solution  explains  the  correlation  between  the  dimensionless  stand-off  distance  and  the  dimen¬ 
sionless  reaction  rate  parameter  previously  observed  by  Hornung  (1972)  for  nitrogen.  The  physics 
of  the  correlation  is  equivalent  to  binary  scaling.  Based  on  the  solution,  a  new  dimensionless 
reaction  rate  parameter  is  defined  to  generalize  Hornung’s  correlation  for  more  complex  gases 
than  nitrogen.  The  physical  meaning  of  fi  is  recognized  as  the  ratio  of  the  energy  absorption  rate 
by  chemical  reactions  just  after  the  normal  shock  to  the  input  rate  of  free-stream  kinetic  energy. 
This  suggests  that  there  exists  a  universal  relationship  which  applies  to  all  gases  and  free-stream 
conditions,  provided  that  the  Mach  number  is  sufficiently  high.  The  form  of  the  function  is  also 
determined  by  the  approximate  analysis  and  is  presented  together  with  computational  results  for 
air  in  Fig.  1.11  and  together  with  experimental  results  from  T5  for  carbon  dioxide  in  Fig.  1.12. 
Experimental  and  numerical  results  confirm  the  new  correlation. 
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Fig.  1.11.  Dimensionless  stand-off  distance  plotted  against  reaction  rate  parameter  for  different  spe¬ 
cific  reservoir  enthalpies  in  air.  The  symbols  represent  computational  results.  The  filled 
circle  is  the  point  at  which  the  curve  branches  off  from  the  solid  line  for  each  pelp»k- 


Q 

Fig.  1.12.  Dimensionless  stand-off  distance  plotted  against  reaction  rate  parameter  for  different  spe¬ 
cific  reservoir  enthalpies  in  carbon  dioxide.  The  symbols  represent  experimental  results. 


The  effect  of  nonequilibrium  recombination  downstream  of  a  curved  two-dimensional  shock 
was  also  investigated.  An  analytical  solution  for  an  ideal  dissociating  gas  was  obtained,  giving 


an  expression  for  dissociation  fraction  as  a  function  of  temperature  on  a  streamline.  The  solution 
agrees  well  with  the  numerical  result  and  provides  a  rule  of  thumb  to  check  the  validity  of  binary 
scaling  for  the  experimental  conditions.  The  effects  upon  the  binary  scaling  of  the  large  difference 
in  free-stream  temperature  between  flight  and  free-piston  shock  tunnel  conditions  are  discovered. 

Though  theoretical  and  computational  results  show  that  recombination  is  important  in  parts 
of  the  flow  field,  binary  scaling  is  a  satisfactory  model  for  global  quantities  such  as  the  stand-off 
distance. 

1.2.3  Vorticity  generation  by  curved  shock 

Fig.  1.13(a)  shows  a  resonantly-enhanced  shadowgraph  of  the  carbon  dioxide  flow.  Stream¬ 
lines  are  marked  by  the  sodium  seeding.  At  the  side  of  the  sphere  these  show  similarities  with 
the  instability  of  a  shear  layer.  Fig.  1.13(b)  shows  these  vortical  structures  on  a  larger  scale.  The 
vortical  structures  may  be  the  the  result  of  Kelvin- Helmholtz  instability  of  the  high-vorticity  layer 
at  the  edge  of  the  high-entropy  layer  near  the  body.  If  they  are,  to  our  knowledge,  this  resonantly- 
enhanced  shadowgraph  is  the  first  visualization  of  this  phenomenon.  Further  experiments  are 
needed  before  definite  conclusions  may  be  reached.  Nevertheless,  the  resonantly-enhanced  shad- 
owgraphic  technique  provides  a  promising  tool  for  studying  the  vorticity-interaction  problem  in 
hypervelocity  blunt-body  flows. 


Fig.  1.13.  Resonantly-enhanced  shadowgraph  of  CO2  flow  over  a  sphere  of  4  in  diameter.  Free- 
stream  condition:  CO2,  3.55  km/s,  0.081  kg/m3,  11.27  MJ/kg. 

The  interaction  of  the  entropy  layer  and  the  boundary  layer  has  been  a  challenging  aerodynamic 
problem  for  years.  The  effect  of  the  shock-generated  free-stream  vorticity  on  hypersonic  boundary 
layers  needs  to  be  understood.  Also,  experimental  and  numerical  investigations  of  the  real-gas 
effects  in  this  important  feature,  both  with  respect  to  the  amount  of  vorticity  produced  by  the 
bow  shock,  and  with  respect  to  the  additional  amount  produced  by  dissociation  downstream  of  the 
bow  shock  are  essential  to  improve  the  understanding  of  the  fundamental  physics  in  nonequilibrium 
flow. 
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CHAPTER  II 

SHOCK- VORTICITY  INTERACTION 

11.1  NONEQUILIBRIUM  CHEMISTRY  IN  SHOCK- VORTEX  INTERACTION 

Objectives  and  Status  of  Research 

The  objectives  of  this  research  are  to  elucidate  the  generic  mechanisms  controlling  (1)  interac¬ 
tion  of  a  shock  wave  with  a  held  of  vortidty  in  the  presence  of  nonequilibrium  chemistry  and  (2) 
vorticity  generation  mechanisms  in  shock  focusing  with  reaction  chemistry.  Progress  has  been  made 
in  computation  of  steady  compressible  vortex  structures  which  will  then  be  used  as  special  initial 
conditions  for  controlled  numerical  experiments  on  shock  vorticity  interaction.  Progress  has  also 
been  made  in  developing  the  computational  capability  for  simulation  of  vorticity  generation  through 
shock  focusing.  A  large  scale  two  dimensional  compressible  flow  solver  has  been  constructed  and  a 
high  resolution  study  has  been  started  on  two  massively  parallel  machines  including  the  Intel  Delta 
and  the  IBM  SP/1.  Finally  we  have  applied  Whitham  theory  to  compute  an  approximation  for 
the  vorticity  generated  by  a  shock  as  it  interacts  with  a  density  gradient.  Computations  are  now 
in  progress  to  test  this  approximation. 

The  next  stage  of  this  work  will  be  the  use  of  the  steady  vortex  configurations  computed 
above  as  initial  conditions  for  our  flow  simulation  and  the  detailed  study  of  the  evolving  vorticity 
distribution.  In  addition  preliminary  explorations  of  the  effect  of  chemistry  will  be  undertaken. 
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(4)  Mr.  Mark  Meloon,  graduate  student.  Thesis  will  entail  numerical  study  of  vorticity  generation 
through  mechanisms  of  shock  focusing,  and  shock-vortex  interactions  with  reaction  chemistry. 

11.1.1  Steady  Compressible  Vortex  Arrays 

As  part  of  our  study  of  shock-vortex  interactions  progress  has  been  made  towards  computing 
steady  compressible  vortex  solutions.  Our  plan  is  to  use  these  solutions  as  initial  conditions  and  to 
then  examine  the  shock-vortex  interaction  resulting  from  the  interaction  of  an  incipient  shock  with 
the  vortex  distribution  both  in  the  presence  and  in  the  absence  of  chemistry  effects.  This  approach 
has  the  advantage  that  the  vortex  distribution  is  steady  and  hence  one  does  not  have  to  take  into 
account  any  time  dependence  prior  to  the  interaction  of  the  shock  with  the  vortex  solution.  At 
present  our  work  has  concentrated  on  two  models  of  steady  compressible  vortices.  The  first  is  an 
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array  of  hollow  core  vortices,  an  extension  of  an  approach  initiated  by  Moore  and  Pullin  ( JFM„ 
185, 188, 1987).  This  type  of  solution  is  that  it  allows  us  to  use  hodograph  techniques  which  make 
it  possible  to  apply  semi-analytic  techniques  to  effect  the  solution. 

This  problem  can  be  formulated  by  use  of  the  stream  function  ip(x,y)  such  that  the  velocity 
field  given  by  u  =  P~l(dip/dy),  v  =  —p~i(drp/dx),  where  p  is  the  local  density.  The  hodograph 
plane  method  then  allows  derivation  of  the  Chaplygin  equation  for  ip 


*  l1 — zr)w+,v — 2T-;97+l1 — w~)w =0- 


where  ca  is  the  stagnation  sound  speed.  We  then  determine  the  image  of  the  row  of  vortices 
in  the  hodograph  plane.  The  main  advantage  of  this  approach  is  that  since  the  velocity  on  the 
vortex  boundaries  must  be  tangent  to  the  boundaries  and  since  the  core  is  stagnant  the  hodograph 
transformation  maps  the  vortex  boundary  to  a  straight  segment. 


We  have  solved  this  problem  using  two  complementary  approaches.  The  first  is  to  use  the 
Rayleigh-Janzen  method  to  generate  a  perturbation  solution  in  powers  of  the  square  of  the  Mach 
number.  This  has  been  completed.  We  have  also  set  up  a  numerical  approach  by  solving  the 
Chaplygin  equation  using  finite  differences.  In  order  to  obtain  high  resolution  we  have  adopted  a 
capacitance  matrix  technique  which  allows  us  to  solve  for  a  solution  which  satisfies  the  boundary 
conditions  described  above  through  the  modification  of  a  solution  which  can  be  obtained  using  a 
fast  Poisson  solver.  At  present  we  are  comparing  the  results  of  this  numerical  approach  along  with 
the  Rayleigh-Janzen  approach. 


The  second  model  is  obtained  from  a  continuation  of  the  analytical  solutions  of  Stuart  ( J.F.M. , 
29, 417, 1967)  for  a  row  of  steady  incompressible  vortices ,  to  finite  Moo.  This  reduces  to  the  solution 
of  coupled  nonlinear  equations  for  the  density  and  stream-functions 


+  PooP2*-2*1'9  =  0, 

-  p-'e-WtP  =  0, 


where  p  =  p/pa,  q  =  q/ca,  the  subscripts  “s”  and  “oo”  refer  to  stagnation  conditions  and  conditions 
at  infinity  respectively,  and  ip  is  a  non-dimensional  stream-function.  With  suitable  boundary  con¬ 
ditions  these  equations  can  be  solved  by  finite  difference  techniques  to  yield  steady  “ compressible 
Stuart  vortex”  solutions  of  the  Euler  equations.  These  can  then  be  used  as  initial  conditions  for 
the  unsteady  shock-vortex  interaction  codes.  This  work  is  in  progress. 


II.1.2  Numerical  simulation  of  shock  generated  vorticity 

Several  codes  for  the  computation  of  2-D  unsteady  Euler  flows  have  been  written,  one  of  which 
( PGP2D:  Parallel  General  Purpose  2D)  is  ported  to  the  Caltech  Intel  Paragon  computer.  We  de¬ 
scribe  numerical  computations  in  two  space  dimensions  of  Richtmyer-Meshkov  instability  initiated 
through  the  interaction  of  a  propagating  shock  and  a  perturbed  density  gradient.  These  calcula¬ 
tions  utilize  the  EFM  numerical  technique  (Pullin  JCP,  34,  231,  1980),  a  truly  multidimensional 


algorithm  for  which  the  addition  of  finite-rate  chemistry  is  straightforward.  The  initial  condition 
is  a  shock  of  given  Mach  number  impinging  on  a  density  gradient  given  by 

representing  a  sinusoidally  perturbed  interface  where  p  is  related  to  the  asymptotic  density  ratio  by 
P  -  (PilPl)  ~  l)/2—  We  take  the  interface  to  be  given  by  zj(p)  =  *o  +  Asin(2irny),  with  A  =  0.1, 
p2/pl  =  2,  n  =  2,  and  6  -  0.01.  To  the  left  of  the  layer  we  place  an  oncoming  shock  with  shock 
Mach  number  of  10.  The  resolution  of  the  calculations  is  400  X  200  points. 

The  results  of  the  calculations  are  displayed  in  the  form  of  color  images  in  Figures  II.l  and 
II.2.  In  Figure  Il.la-d  we  display  a  time  sequence  of  the  density.  The  calculations  clearly  show  the 
diffraction  of  the  shock  as  it  interacts  with  the  heavier  medium  down  stream.  A  reflected  shock  is 
formed  as  well  which  itself  focuses  as  the  flow  develops.  The  density  gradient  displays  Richtmyer- 
Meshkov  instability  and  growth  is  observed  as  the  layer  propagates  downstream.  Figure  II.2a-d 
displays  the  corresponding  vorticity.  One  very  interesting  feature  is  the  production  of  vortidty  of 
opposite  sign  in  the  later  stages  of  the  calculation.  This  is  clearly  seen  for  example  in  Figure  II.2d 
where  vorticity  of  opposite  sign  to  the  primary  vorticity  deposited  by  the  initial  interaction  of  the 
shock  with  the  density  gradient  is  formed  ahead  of  the  density  gradient.  This  is  presumably  due 
to  shock  oscillations  of  the  type  seen  in  our  calculations  of  shock  dynamics. 

We  have  also  derived  a  formula  to  approximate  the  vorticity  generated  by  the  shock  using 
Whitham’s  shock  dynamics.  Associated  with  the  formulation  of  geometrical  shock  dynamics  is  the 
concept  of  a  local  flow  behind  a  shock  element.  We  consider  the  flow  at  two  closely  spaced  times  t 
and  t  +  At.  The  idea  is  to  use  the  formula  relating  circulation  to  vorticity  as  follows: 


where  u  is  the  vorticity  contained  within  the  cell,  A  is  the  area  of  the  cell,  u  is  the  local  velocity 
along  the  boundaries  of  the  cell,  and  1  represents  the  element  of  arc  along  the  cell  boundary.  In  using 
the  formula  above  we  assume  that  the  velocities  at  times  t  +  At  have  not  changed  appreciably  from 
their  previous  values.  This  is  strictly  speaking  not  true  and  in  order  to  improve  this  a  better  model 
of  the  flow  behind  the  shock  needs  to  be  developed.  The  theory  of  geometrical  shock  dynamics 
considers  only  local  fluid  effects  in  order  to  determine  the  motion  of  the  shock.  An  improved  theory 
has  been  developed  recently  by  J.  P.  Best  and  it  would  be  of  interest  to  determine  if  this  could 
be  used  to  provide  improved  values  of  the  vorticity,  but  this  has  not  been  implemented  here.  The 
result  of  this  calculation  is 

u_  1 

aoM  de  y  7  +  1  J 

where  ao  is  the  local  sound  speed  ahead  of  the  shock  and  M  is  the  local  shock  Mach  number.  Note 
that  this  result  implies  there  are  two  ways  in  which  vorticity  is  generated.  The  first  is  through 
the  variation  of  the  sound  speed  along  the  shock.  This  would  occur  for  example  if  the  shock 
encounters  a  density  gradient.  The  second  is  through  variations  of  the  Mach  number  along  the 
shock.  This  implies  that  the  shock  is  curved.  Both  these  mechanisms  are  known  to  generate  vorticity 
in  compressible  flow  through  the  interaction  of  the  misaligned  density  and  pressure  gradients.  Note 
also  that  no  vorticity  is  generated  in  the  absence  of  these  conditions.  Further  computations  are 
now  in  progress  to  compare  the  actual  vorticity  levels  predicted  by  shock  dynamics  with  those  of 
the  full  calculation. 
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(a)  t  =  0.03 


(b)  t  =  0.06 


(c)  t  =  0.09 


(d)  t  =  0.12 
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Figure  II. 1:  Shock-density  gradient  interaction,  M=10. 
Time  sequence  of  color  coded  density  field 
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(d)  t  =  0.12 


Figure  II.2:  Shock-density  gradient  interaction,  M=10. 
Time  sequence  of  color  coded  vorticity  field 


II.2  SHOCK  WAVE  INTERACTIONS  IN  HYPERVELOCITY  PLOW 


Objectives  and  Status  of  Research 

The  objective  of  this  experimental  study  is  to  determine  the  effects  of  chemical  relaxation 
on  the  shock-on-shock  problem,  and,  in  particular,  to  quantitatively  demonstrate  the  magnitude 
of  nonequilibrium  effects  on  shock  impingement  heating.  Preliminary  experiments  in  T5  were 
described  in  the  proposal  for  this  grant  and  in  Ref.  1.  Presently,  cylindrical  and  spherical  models 
are  being  constructed  for  a  new  series  of  tests.  Modifications  are  required  to  the  impinging-shock 
generators  (items  B.l  and  2,  Statement  of  Work;  Design  and  construct  shock  generators),  and  the 
three  T5  test  conditions  established  for  this  project  are  being  refined. 

The  primary  effort  during  this  reporting  period  was  to  make  major  improvements  in  methods 
for  heat  transfer  measurement  and  optical  interferometry  in  T5.  In  addressing  item  A.3,  Statement 
of  Work  (Improve  existing  Wollaston  prism  differential  interferometer),  a  more  sweeping  solution 
was  sought;  a  completely  new  holographic  system  for  optical  diagnostics  in  T5  has  been  designed 
and  built.  This  instrument  will  be  reported  at  another  time.  In  this  report  we  document  the 
development  of  a  new  heat  transfer  gauge  which  gives  performance  superior  to  any  existing  device. 


II.2.1  A  surface  junction  thermocouple  sensor  for  transient  heat  flux  measurements 

The  measurement  of  heat  transfer  rates  to  aerodynamic  surfaces  is  important  in  the  design  of 
thermal  protection  for  hypervelocity  aerospace  vehicles.  These  flows  may  be  created  in  wind  tunnels 
for  only  a  fraction  of  a  second,  due  to  the  high  specific  energies  involved  and  the  limited  strength 
of  materials  at  elevated  temperatures.  Thus  the  problem  is  transient  in  nature  and  a  common 
measurement  technique  is  to  sense  the  time  history  of  the  surface  temperature  of  the  model.  The 
unsteady  conduction  of  heat  within  a  solid  body  is  described  by  linear  partial  differential  equations 
and  so  the  heat  transfer  rates  at  the  surface  of  a  body  may  be  inferred  from  the  surface  temperature 
measurements.  The  operating  principles,  and  typical  sensor  designs  are  described  by  Schultz  and 
Jones  (1973). 

We  have  conducted  experiments  in  the  T5  free  piston  driver  shock  tunnel  at  GALCIT.  This 
facility  is  capable  of  simulating  the  high  density  hypervelocity  flows  about  space  vehicles  at  speeds 
of  order  5000  m/s.  Under  these  extreme  conditions  the  flow  in  typical  facilities  is  contaminated 
with  the  products  of  mechanical,  physical  and  chemical  breakdown  of  a  variety  of  metallic  and 
nonmetallic  materials.  Some  erosion  of  the  surface  of  the  gauge  is  probable  under  these  conditions 
and  this  impacts  the  design  of  the  transducer. 

A  commonly  used  technique  for  gauge  fabrication  is  to  deposit  a  thin  metal  film  (typically 
Platinum)  on  a  glass  or  ceramic  substrate.  The  surface  temperature  is  then  sensed  through  the 
temperature  dependence  of  the  resistivity  of  the  metal  film.  This  technique  yields  a  relatively 
sensitive  device  however  the  thin  film  is  prone  to  damage  in  hostile  environments.  Any  erosion  of 
the  metal  film  alters  the  sensitivity  of  the  gauge  and  so  these  devices  are  unsuitable  for  some  wind 
tunnels.  Observe  also  that  since  the  sensitivity  of  the  gauge  depends  on  the  geometry  of  the  device, 
calibration  of  each  individual  sensor  is  required. 

An  alternative  physical  mechanism  for  temperature  measurement  is  the  thermoelectric  EMF 
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produced  at  a  junction  of  dissimilar  metals.  If  a  metallic  substrate  is  used,  the  gauge  is  capable  of 
tolerating  a  larger  transient  heat  flux  due  to  the  higher  thermal  diffusivity  and  the  larger  allowable 
temperature.  The  figure  of  merit  is  TmaxVpck  Eqn.  (4)].  A  particularly  robust  design  results 
when  the  thermocouple  material  itself  is  used  as  the  substrate.  The  sensitivity  of  the  gauge  is  then 
unaffected  by  surface  erosion  or  gauge  geometry  and  depends  only  on  physical  properties  of  the 
materials.  Thus  calibration  of  individual  gauges  is  not  required. 

The  design  problem  in  producing  such  a  gauge  is  to  insulate  the  two  conductors  that  form 
the  thermocouple  so  that  only  a  thin  surface  junction  is  formed.  Common  practice  is  to  construct 
the  gauge  from  coaxial  cylindrical  thermocouple  elements  that  are  insulated  by  a  layer  of  epoxy 
or  ceramic  adhesive.  The  surface  junction  is  formed  by  electrolytic  deposition  of  a  metal  film  or 
by  abrasion  of  the  exposed  surface.  In  either  case  very  fine  tolerances  are  required  to  produce  a 
small  gap  between  the  thermocouple  elements  and  the  effective  depth  of  the  junction  scales  with 
the  width  of  the  gap  between  the  conductors.  Note  that  as  the  test  time  is  reduced,  a  thinner 
junction  is  required  for  the  same  materials. 


II.2.2  Analysis  of  heat  conduction  in  gauge  substrate 

In  order  to  illustrate  the  principles  of  gauge  operation  and  to  assess  the  effect  of  the  finite  thickness 
of  the  surface  junction,  consider  the  unsteady  conduction  of  heat  in  a  one-dimensional  semi-infinite 
solid.  For  small  temperature  gradients  this  process  is  described  by  the  usual  linear  partial  differ¬ 
ential  equation, 

d2T  _  ldT  _  k 

dx 2  a  dt'  °  pc'  * 

where  T  is  the  temperature,  x  is  the  spatial  coordinate  normal  to  the  surface,  t  is  time  and  a  is 
the  thermal  diffusivity,  k  is  the  thermal  conductivity,  p  is  the  density  and  c  is  the  specific  heat  of 
the  material.  For  the  following  uniform  initial  condition,  7i,  and  instantaneously  applied,  constant 
surface  heat  flux,  qo; 

T(x,0)  =  Th 


(dT\  _-9o 

V0*;*=o-  *  ’ 


t  >  0, 


the  solution  is, 


Expanding  this  solution  asymptotically  in  powers  of  t  as  t  — *  oo  at  constant  x  we  obtain. 

Further,  the  first  term  is  asymptotically  dominant  for  *  <  <<  1. 

(«*)* 

k  y  x 
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and  this  expresses  the  usual  result  that  the  surface  temperature  rises  parabolically  with  time  for  a 
constant  heat  flux.  The  criterion  — <<  1  provides  a  limit  on  the  thickness  of  the  surface  sensing 

M* 

element  for  practical  designs.  The  next  order  term,  ^x,  gives  the  error  in  the  surface  temperature 
measurement  associated  with  this  finite  thickness  and  the  estimate  is  also  asymptotically  valid 
provided  — *-r  <<  1.  Extracting  the  exponential  prefactor  and  expanding  asymptotically  as  t  — *  0 

(ari)J 

at  constant  x,  the  leading  order  behavior  is; 


(5) 


Hence  at  small  times  the  temperature  rises  exponentially  fast.  The  significance  of  the  various 
terms  of  the  expansion  is  illustrated  in  Fig.  II.3  At  long  times  the  solution  asymptotes  to  the  true 
surface  temperature,  given  by  the  first  term  in  the  expansion  (3),  with  a  constant  error  given  by 
the  second  term  of  (3).  This  error  becomes  small  with  respect  to  the  true  surface  temperature 
provided  -*-r  <<  1. 

(•«)* 


Fig.  II.3.  —  Asymptotic  behavior  of  a  surface  temperature  sensor.  For  illustrative  purposes  a  =  1, 
fc  =  1,  90  =  1,  x  =  1. 


II. 2. 3  Description  of  new  gauge  design 

A  new  gauge  design  has  been  developed  and  tested  with  the  aim  of  simplifying  the  construction  of 
surface  junction  thermocouple  sensors.  The  essential  feature  of  this  design  is  the  use  of  a  slightly 
tapered  center  conductor.  The  resulting  interference  between  the  tapered  center  pin  and  the  sharp 
edged  outer  conductor  forms  a  thin  surface  junction  and  provides  reliable  insulation  (see  Fig.  II.4). 
The  protruding  pin  is  severed  and  the  surface  is  polished  once  the  gauge  has  been  assembled. 


Fig.  II.4.  Layout  of  new  surface  junction  thermocouple  sensor. 


This  design  results  in  a  robust  transducer  whose  sensitivity  is  determined  solely  by  the  physical 
properties  of  the  thermocouple  materials.  There  is  a  further  weak  dependence  on  the  taper  angle 
however  this  is  constant  for  all  gauges  and  is  determinable  by  either  calibration  or  analysis. 

The  effective  thickness  of  the  junction  may  be  estimated  from  the  surface  roughness  of  the 
thermocouple  components  and  the  taper  angle  of  the  center  conductor  (see  Fig.  II.5).  The  thickness 
of  the  junction,  averaged  over  the  two  surfaces  may  be  estimated  to  be, 

thickness  =  taper  x  roughness. 

It  is  possible  to  achieve  a  0.3  jinish  with  conventional  machining  operations  and  with  a  5:1 
taper  the  average  junction  depth  is  1.5  /im.  For  the  gauges  described  here  type  E  (Chromel  - 
Cons  tan  tan)  thermocouples  were  used  and  the  average  thermal  diffusivity  for  these  materials  is 
5  X  10-6  m2/a.  A  response  time  of  0.5  us  then  follows  from  the  condition  — *-»-  <<  1.  The  test 

(at)* 

time  for  typical  hypervelocity  wind  tunnels  lies  in  the  range  1-100  ms  and  so  this  is  acceptable. 
The  length  of  the  gauge  must  be  sufficient  such  that  the  temperature  rise  at  the  rear  surface  is 

negligible.  This  is  determined  from  the  condition  — *-r  »  1;  where  here  the  scales  are  understood 

(at)* 

to  be  the  test  time  and  the  length  of  the  gauge. 
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A  second  dissimilar  metal  junction  is  formed  between  the  outer  conductor  and  the  model. 
This  junction  extends  along  the  full  length  of  the  gauge  since  these  surfaces  are  not  insulated  and 
so  the  thermoelectric  EMF  is  small.  This  common  mode  EMF  and  noise  pickup  are  reduced  by 
the  use  of  a  relatively  crude  differential  amplifier  (see  Fig.  II.6).  The  cold  junction  is  formed  at 
ambient  temperature  at  the  circuit  board.  This  arrangement  is  adequate  given  the  constancy  of 
the  cold  junction  temperature  during  the  short  run  time  and  the  large  hot  junction  temperature 
rise  experienced  in  hypervelocity  wind  tunnels  (approx.  200K). 


Fig.  11.6.  Schematic  of  amplifier  circuit. 


The  center  pin  is  located  in  the  concentric  outer  conductor  and  secured  with  epoxy  or  ceramic 
adhesive.  Sufficient  force  is  applied  to  ensure  adequate  registration  of  the  mating  parts  however 
care  must  be  taken  not  to  deform  the  edge  of  the  outer  conductor.  Thermocouple  extension  wires 
are  resistance  welded  to  the  rear  of  the  assembly  by  discharging  a  capacitor  through  the  wire. 
The  thermocouple  is  then  installed  flush  with  the  outer  surface  of  the  model  and  the  protruding 
center  pin  is  removed.  The  excess  material  is  removed  by  filing  and  the  surface  is  honed  to  produce 
the  surface  junction.  It  is  helpful  to  apply  a  small  amount  of  Prussian  blue  to  the  surface  of  the 
thermocouple  and  once  this  is  just  removed  by  filing,  the  formation  of  a  thin  junction  is  assured. 
It  is  advisable  to  manufacture  the  model  from  a  material  somewhat  harder  than  the  thermocouple 
material  so  that  the  model  contours  are  not  altered  during  this  process.  It  is  possible  for  test  section 
models  to  experience  surface  temperature  rises  of  1000K  during  the  short  test  time.  Degradation 
of  the  adhesive  bond  between  the  conductors  will  also  occur  however  the  penetration  depth  during 
the  test  time  is  small  relative  to  the  overall  length  of  the  gauge  and  the  bond  at  the  rear  maintains 
the  integrity  of  the  device. 

II.2.4  Deducing  the  surface  heat  transfer  rate 

If  the  unsteady  conduction  of  heat  in  the  thermocouple  sensor  is  modeled  by  the  analysis  presented 
in  §11.2.2  then  given  the  time  history  of  the  surface  temperature  and  the  physical  properties  of  the 
gauge  materials  it  is  possible  to  solve  an  inverse  problem  for  the  time  history  of  the  surface  heat 
flux.  The  procedure  is  described  by  Schultz  and  Jones  (1973)  whereby  the  calculation  is  performed 
using  a  problem  specific  quadrature  scheme  or  by  an  electronic  analogue. 

Observe  however  that  the  governing  equations  are  linear  and  hence  the  dynamic  behavior  of 
such  a  system  may  be  generally  be  described  by  a  convolution  integral, 

x,t)=  f  g(x,t  —  r)g(r)dr. 

JO 


(6) 
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The  unit  impulse  response  function,  g(z,  t),  is  given  by  the  derivative  of  Eqn.  (2), 


?(*»  0 


dAT(x,t) 

3t 


t  >  0. 


(7) 


Note  that  the  integral  (6)  is  improper  in  the  limit  x  -*  0  and  a  finite  value  of  x  is  taken  to  regularise 
the  problem  (as  estimated  in  §11.2.2).  Recasting  the  problem  in  this  general  form  allows  the  use  of 
spectral  deconvolution  methods  that  provide  a  mean  square  optimal  estimate  of  the  time  history 
of  the  surface  heat  flux.  This  procedure  provides  optimal  rejection  of  the  errors  associated  with 
the  differentiation  of  noisy  experimental  data.  Theoretical  background  is  given  in  Papoulis  (1984) 
and  the  numerical  implementation  with  discrete  Fourier  transform  algorithms  is  described  by  Press 
et  al  (1988).  A  concise  description  of  these  methods  applied  to  an  analogous  problem  is  given  by 
Sanderson  and  Simmons  (1991). 


If  the  noisy  experimental  data  in  discrete  time  is. 


Zj  =  A  Ti  +  n„ 


(7) 


then  the  deconvolution  problem  comprises  determining  the  acausal  filter,  <f>i,  that  gives  a  mean 
square  optimal  estimate,  q* ,  of  the  time  history  of  the  surface  heat  flux.  Thus, 

N- 1  . 

5^  4>i-kZk  =  9i-kik ;  « =  o»  - 1.  (8) 

k—0  k= 0 

By  the  convolution  property  of  the  discrete  Fourier  transform,  and  since  g,  =  0  for  i  <  0, 


SjZj^GjQ)  j  =  0, 1...JV  —  1. 


It  may  be  shown  that  the  optimal  filter  is  approximated  by, 


SzZj  ~  vj 

$ZZj 


(9) 


(10) 


where  SzZj  >s  the  discrete  power  spectral  density  of  the  surface  temperature  measurement  and  Vj 
is  the  power  spectral  density  of  the  signal  noise.  By  inverse  discrete  Fourier  transform,  the  optimal 
estimate  of  the  surface  heat  flux  is  then, 


7*  =  F"1 


9. 


i,j  =  0,1. ..N  -  1. 


(11) 


Note  that  all  time  series  must  be  padded  to  twice  their  length  with  trailing  zeroes  to  enforce  causality 
of  the  right  hand  side  of  Eqn.  (8)  when  discrete  Fourier  transforms  are  used  in  this  manner. 


II.2.5  Experimental  results 

In  order  to  demonstrate  the  performance  of  the  new  heat  transfer  gauge  we  have  measured  the 
stagnation  point  heat  transfer  rate  to  a  50 mm  diameter  sphere  in  hypervelocity  flow.  Fig.  II.7 
shows  the  time  history  of  the  nozzle  reservoir  pressure,  test  section  pitot  pressure,  stagnation  point 
surface  temperature  and  the  stagnation  point  heat  transfer  rate  deduced  from  it. 
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Fig.  II.7.  Stagnation  point  heat  transfer  measurement  for  sphere  in  hypervelocity  flow.  The  free 
stream  conditions  are  p  =  0.020kg/m3,  u  =  4590m/s,  T  =  1700K,  p  =  10.4kPa,  M  =  5.65, 
ho  =  13.2MJ/kg  and  these  are  computed  from  measurements  of  the  shock  tube  Ailing 
pressure,  primary  shock  speed  before  reflection  and  the  nozzle  reservoir  pressure. 

11.2.5-1  Accuracy  of  measurements  —  general  considerations 

The  repeatability  of  the  measurement  has  been  established  by  repeating  the  condition  correspond¬ 
ing  to  Fig.  II. 7  with  a  different  gauge  of  the  same  design  and  these  results  are  illustrated  in  Fig.  II.8. 
The  flow  environment  is  extremely  hostile  and  since  the  thermocouple  junction  is  electrically  ex¬ 
posed  to  the  flow  it  is  necessary  to  demonstrate  that  erroneous  EMFb  are  not  generated  by  physical 
mechanisms  other  than  the  thermoelectric  effect.  A  dummy  gauge  of  identical  design  was  man¬ 
ufactured  with  both  conductors  made  of  Constantan  and  so  no  thermoelectric  EMF  should  be 
generated.  Fig.  II.9  shows  the  results  obtained  with  this  dummy  gauge  at  a  high  enthalpy  condi¬ 
tion.  Apart  from  an  initial  transient  that  decays  rapidly,  no  erroneous  EMFs  are  discernible  on  this 
scale.  The  deconvolved  heat  transfer  rate  is  three  orders  of  magnitude  less  than  the  actual  heat 
transfer  rate. 

The  conversion  from  thermoelectric  EMF  to  temperature  is  given  by  the  standards  for  type- 
E  thermocouples  wl.hin  prescribed  error  limits.  Further  errors  associated  with  the  uncertainty 
in  the  physical  properties  of  the  thermocouple  elements,  the  difference  in  a  and  k  for  the  two 
dissimilar  metal  conductors  and  the  reduced  heat  conduction  away  from  the  surface  due  to  tapering 
of  the  center  conductor  may  be  eliminated  by  calibration.  The  tapering  of  the  Constantan  pin 
effectively  reduces  its  thermal  diffusivity  and  this  error  tends  to  negate  the  unavoidably  lower 
thermal  diffusivity  of  Chromel.  These  systematic  errors  may  be  eliminated  by  calibration  to  any 
desired  level  of  accuracy.  The  remaining  systematic  error  due  to  the  finite  thickness  of  the  surface 
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Fig.  II  8.  - ' .  ■*  -\tion  point  heat  transfer  measurement  for  sphere  in  hypervelocity  flow.  The  free 
.»»  <*.  jn  conditions  are  nominally  identical  with  those  for  Fig.  II.8 

junction  may  be  estimated  from  the  preceding  analysis. 

11.2.5.2  Accuracy  of  measurements  —  high  heat  transfer  rates 

The  surface  temperature  rise  of  models  in  T5  may  easily  exceed  several  hundred  Kelvin.  Under 
these  conditions  the  physical  properties  of  the  thermocouple  materials  change  significantly  with 
temperature.  Particularly  strong  variations  of  the  specific  heat  are  possible  in  the  vicinity  of 
solid  state  phase  transitions.  Metastable  phases  that  may  arise  under  conditions  of  rapid  heating 
and  cooling  pose  further  difficulties.  Existing  analyses  for  deducing  the  model  heat  transfer  rate 
from  the  surface  temperature  time  history  are  based  on  linear  hypotheses.  Some  literature  exists 
for  estimating  the  errors  associated  with  the  quasi-linear  temperature  dependence  of  the  material 
properties  [Schultz  and  Jones  (1973),  Hartunian  and  Varwig  (1962)].  The  issue  of  genuine  non¬ 
linearity  associated  with  phase  transitions  remains  open.  The  lack  of  reliable  data  and  the  difficulty 
of  obtaining  it  through  calibration  pose  further  problems.  Consequently,  analysis  of  these  effects 
will  initially  be  limited  to  an  estimation  of  errors. 

II.2.6  Conclusions 

A  new  surface  junction  heat  transfer  gauge  has  been  developed  and  its  use  demonstrated  by  the 
measurement  of  the  stagnation  point  heat  flux  for  a  sphere  in  a  reacting  supersonic  flow.  The 
repeatability  and  validity  of  the  measurements  in  a  hostile  flow  environment  has  been  demonstrated. 
Analysis  has  illustrated  the  principle  of  operation  and  spectral  methods  were  used  to  infer  the 
surface  heat  flux  from  surface  temperature  measurements.  Potential  error  sources  were  identified. 


Tim«  (mi)  Tim«  (mi) 

Fig.  II.9.  Stagnation  point  heat  transfer  measurement  for  sphere  in  hypervelocity  flow,  using  dummv 
Constantan  -  Constantan  thermocouple.  The  free  stream  conditions  are  p  =  0.018kg/m  , 
u  =  5300m/s,  T  =  2320K,  p  —  13.7kPa,  M  =  5.52,  ho  =  19.6MJ/kg  and  these  are 
computed  from  measurements  of  the  shock  tube  filling  pressure,  primary  shock  speed 
before  reflection  and  the  nozzle  reservoir  pressure.  (The  large  drop  in  the  pitot  pressure 
at  t  =  1.5ms  is  intentionally  caused  by  the  passage  of  an  oblique  shock  wave.) 
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CHAPTER  III 

SUPERSONIC  SHEAR-LAYER  FLOWS 


III.I  Bi-supersonic  shear-layer  flows 

During  the  past  year,  the  GALCIT  Supersonic  Shear  Layer  (S3L)  facility  was  upgraded  to 
permit  operation  in  a  bi-supersonic  flow  regime,  i.e.,  both  freestreams  supersonic.  Investigation  of 
this  fully-hyperbolic  flow  regime  will  allow  us  to  address  a  host  of  important  issues,  such  as  the 
influence  of  upstream  and  downstream  boundary  conditions,  the  interaction  of  shock  waves  with  the 
turbulent  mixing  zone,  and  the  effect  of  new  compressibility  phenomena  in  both  chemically-reacting 
and  non-reacting  gas-phase  flows. 


r 


Fig.  HI.l.  Subsonic/supersonic  shear  layer  schlieren-image  data.  Inlet  conditions:  M\  =  0.65  [N2], 
A/2  =  1.13  [N2] 

The  first  flow  condition  investigated  in  the  upgraded  facility,  a  M\  =  0.65  [N2J,  M2  =  1.13  [N2] 
(inlet  conditions)  shear  layer,  is  shown  in  the  form  of  a  schlieren  photograph  visualization  in  III.I. 
Although  this  flow  is  not  bi-supersonic  at  the  beginning  of  the  test  section  (left  side  of  the  im¬ 
age),  a  fan  produced  by  a  5°  expansion  turn  in  the  lower  freestream,  on  the  bottom  wall,  just 
upstream  of  the  inlet,  curves  the  shear  layer.  This  results  in  an  aerodynamic  Laval  nozzle  throat, 
in  a  manner  which  causes  a  Mach  number  increase  in  both  streams.  This  produces  bi-supersonic 
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flow  downstream  of  the  aerodynamic  throat,  providing  an  interesting  combination  of  mixed  (liv- 
perbolir/elliptic)  boundary  conditions  for  this  flow  configuration. 

The  first  bi-supersonic  condition  investigated  was  of  a  shear-layer  flow  with  inlet  conditions 
M\  =  1.50  [No],  A/o  =  1.13  [No].  A  schlieren  photograph  appears  in  III. 2.  The  most  conspicuous 
difference  between  this  flow,  as  well  as  the  downst  ream  portion  in  II 1. 1 .  and  the  supersonic/subsonic 
flows  that  had  been  previously  investigated  in  this  facility,  i.f.,  supersonic  high-speed  stream  and 
subsonic  low-speed  stream  (Hall  ft  at.  1993),  is  the  transmission  of  the  high-speed  stream  oblique 
shock  waves  into  the  low-speed  st  ream.  Such  a  t  ransmission  of  stationary  (in  the  laboratory  frame) 
shocks  through  the  turbulent  mixing  zone  can  only  occur  in  a  bi-supersonic  configuration.  This 
configuration  admits  rather  complex  interactions,  as  evidenced  by  the  four-wave  system  observed 
just  above  the  rightmost  numeral  "4"  in  the  run  identification  number,  “404”. 

;,1  rom  the  schlieren  photograph,  it.  can  be  seen  that  the  visual  grov  '■  rate  of  the  shear  layer 
is  not  substantially  altered  by  these  wave-turbulence  interactions.  This  insensitivity,  however,  is 
to  be  expected,  given  the  weak  nature  of  the  incident  waves  and  the  near-unity  density  ratio 
(P'i/P I  ~  0.89)  of  this  flow.  Additionally,  the  shear  layer  growth  rate  realized  in  this  experiment 
is  very  nearly  that  of  the  corresponding,  i.f.,  same  velocity  and  density  ratio,  incompressible  shear 
layer.  In  terms  of  the  total  ((invective  Mach  number, 

L'\  -  i ’2 


Me  = 


(8) 


«|  +  «2 

( (•(/..  Sandham  &  Reynolds  1989).  this  result  is  not.  surprising;  this  shear  layer  should  be  regarded 
as  incompressible,  with  A/r  %  0.1  . 


An  additional  feature  of  the  wave-turbulence  interactions  observed  in  this  flow  is  that  shocks 
refract  across  the  shear  layer  almost  discontinuously.  The  waves  exhibit  ukinks”  within  the  shear 
layer.  This  indicates  that  the  transverse  (vertical  in  III.2)  extent  of  the  Mach  number  fluctuations 
is  much  smaller  than  that  of  the  density  fluctuations,  implying  considerably  different  transverse 
extents  (widths)  for  the  Mach  number  and  density  profiles.  Such  behavior  was  also  documented  in 
the  incompressible  shear-layer  investigations  by  Brown  &  Roshko  (1974).  This  difference  may  be  of 
considerable  importance  in  viewing  the  results  obtained  by  linear  stability  analyses,  for  example. 
The  majority  of  such  analyses  employ  Mach  number  and  density  profiles  of  the  same  width  as 
inputs  to  the  computations. 

This  work  is  part  of  the  Ph.D.  research  effort  of  M.  Slessor. 


III. 2  Planar  laser  Rayleigh-scattering  imaging 

As  part  of  our  effort  to  investigate  compressible,  turbulent,  shear-layer  flows  we  have  continued 
the  development  of  planar,  laser-imaging  techniques.  In  recent  experiments,  we  have  extended 
our  previous  measurements  (Rosemann  et  al.  1992)  with  improvements  in  image  signal-to-noise 
ratio.  New  time-resolved  (r  ~  9  ns),  planar,  laser- Rayleigh  scattering  images  have  been  recorded, 
in  non-reacting  supersonic  shear  layers,  for  convective  Mach  numbers  with  respect  to  the  low-speed 
stream,  i.e., 
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in  the  range  of  0.5  <  Mc2  <  1.5.  The  more  recent  results  were  presented  at  the  1993  Annual 
Meeting  of  the  American  Physical  Society,  Division  of  Fluid  Dynamics  (Fourguette  el  al.  1993). 

These  experiments  were  conducted  in  the  Supersonic  Combustion  Facility  in  flow  conditions 
that  had  previously  been  documented.  The  beam  of  a  frequency-doubled  Nd:YAG  laser  (250  mJ 
in  a  9  ns  pulse)  was  steered  through  an  optical  window  into  the  experimental  test-section  and 
shaped  into  a  sheet  of  light.  The  sheet  was  passed  through  slits  in  the  guide  walls  that  acted  as 
“aerodynamic”  optical  windows  and  intersected  the  shear  layer  such  that  a  vector  contained  in  the 
image  plane  was  colinear  with  the  flow  direction.  The  Rayleigh-scattered  light  was  collected  at 
right  angles  to  the  planar  probe  region  and  imaged  onto  a  cryogenically-cooled,  1024  x  1024  pixel, 
CCD  array  of  a  PHOTOMETRICS  camera  system.  The  image  data  were  digitized  using  a  14-bit 
A/D  converter  and  transferred  to  a  VAX  cluster  for  subsequent  storage  and  processing. 

Figures  III.3  and  III.4  depict  examples  of  such  planar-image  data.  The  flows  correspond  to 
shear  layers  with  total  convective  Mach  numbers  (c/.  8)  of  Mc  =  0.54  and  Mc  =  0.94,  respectively. 
In  these  data,  the  flow  is  from  left  to  right.  The  high-speed  stream  occupies  the  upper  part  of  the 
image  and  the  low-speed  stream  the  lower  part  of  the  image.  The  imaged  region  corresponds  to 
a  3.2  X  3.2  cm2  area  in  the  flow,  i.e.,  roughly  twice  the  width  of  the  shear  layer  at  that  location 
(17  cm  downstream  from  the  tip  of  the  splitter  plate). 

In  III.3,  the  high-speed  stream  fluid  is  nitrogen.  In  III.4,  the  high-speed  side  stream  fluid  is 
helium.  In  both  cases,  the  low-speed  stream  fluid  is  ethylene.  The  ratio  of  optical  cross-sections 
(low-speed  fluid/high-speed  fluid)  is  3.79  and  234,  for  the  two  images,  respectively. 

These  data  provide  direct  confirmation  that  compressibility  plays  an  important  role  in  the 
structure  of  the  shear  layer,  in  general,  and  the  mixing  process  in  particular.  Even  at  relatively 
low  compressibility,  i.e.,  Mc  =  0.54  (III.3),  the  mixing  layer  exhibits  a  different  structure  than 


Fig.  III. 3.  Supersonic-subsonic  case:  Mc  =  0.54. 


the  large-scale  roller-like  structures  observed  in  incompressible  shear  layers.  The  low-speed  fluid  is 
seen  to  be  peeled  off  by  the  high-speed  stream.  The  transport  of  fluid  across  the  layer  appears  to 
occur  more  through  relatively  well-defined  intrusions  of  fluid  into  the  opposite  stream,  rather  than 
through  an  entrainment  process  through  the  rotational  motion  of  organized,  large-scale  vortical 
structures.  Nevertheless,  both  these  two  possibilities  are  encountered,  however,  as  reflected  in 
local,  instantaneous  scalar  concentration  profiles,  derived  from  such  images.  These  show  both 
gradual  transitions  from  high-speed  stream  fluid  to  low-speed  stream  fluid,  as  was  documented  by 
Clemens  et  al.  (1991).  We  also  find  "plateaus”,  however,  corresponding  to  regions  of  well-mixed 
fluid  at  an  intermediate  ratio  of  high-  to  low-speed  fluid,  that  are  reminiscent  of  incompressible 
shear-layer  large-scale  vortical-structure  cores. 

In  the  higher  convective  Mach  number  image  data  in  III.4,  the  structure  of  the  shear  layer 
is  different  yet  from  the  previous,  lower  compressibility,  case.  At  these  Mach  numbers,  the  low- 
speed  side  appears  to  be  shredded  by  the  high-speed  stream  and  regions  of  uniformly-mixed  fluid, 
encountered  in  the  Mc  =  0.54  flow,  are  now  not  observed.  In  addition,  a  wave  system  can  be 
clearly  observed  in  the  low-speed  stream.  These  waves  were  observed  in  the  same  facility  using 
Schlieren  imaging  (Hall  et  al.  1991.  Hall  1991).  The  signal-to-noise  ratio  and  resolution  in  these 
images  permits  the  direct  recording  of  these  waves  and  allows  us  to  conclude  that  it  is  the  large- 
scale  structures  that  are  responsible  for  these  waves.  In  particular,  this  wave  system  can  be  seen  to 
be  generated  by  large-scale  structures  in  the  shear  layer  that  are  convecting  with  a  velocity  higher 
that  the  speed  of  sound  in  the  low-speed  stream,  i.e..  an  MC2  >  1-  A  directly-estimated  convective 


Fig.  III.4.  Supersonic-subsonic  case:  Mc  =  0.94. 


Mach  number,  with  respect  to  the  low-speed  stream,  of  A/co  ~  2.5  (9)  for  this  wave  can  be  inferred 
from  the  wave  angles  (/x  ~  24°)  in  the  low-speed  free  stream.  This  estimate  compares  well  with 
the  Hall  et  al.  (1991)  findings  which  were  based  on  schlieren  data  as  well  as  cross-correlation  data 
derived  from  wall-mounted,  high-speed  pressure  transducers. 

This  part  of  the  effort  is  cosponsored  by  the  AFOSR  Air-Breathing  Propulsion  Program.  Con¬ 
tract  No.  F49620-92-J-0290. 

III. 3  Computation  of  non-reacting  flows 

The  new  multidimensional  Godunov  scheme  (Lappas  1993.  Lappas  et  al.  1994a.  1994b)  was 
used  to  compute  a  two-dimensional,  inviscid.  compressible  shear  layer.  The  scheme  employs  the 
method  of  Riemann  Invariant  Manifolds,  which  can  be  viewed  as  the  extension  of  the  method  of 
characteristics  to  the  general  case  of  unsteady,  multidimensional  flow.  Computations  were  per¬ 
formed  to  compare  with  recent,  high-subsonic/supersonic  shear-layer  experimental  results  obtained 
from  the  GALCIT  supersonic  shear  layer  (S3L)  facility.  The  parameters  for  the  computations  are 
chosen  so  that  a  direct  comparison  with  experimental  results  is  possible. 

In  all  the  calculations,  a  perfect  gas  is  assumed,  with  constant,  specific-heat  ratio  of  '■  -  1.4. 
The  computational  domain  is  the  rectangular  region:  0  <  x  <  2  and  0  <  y  <  1.  The  lower 
boundary,  y  —  0.  and  upper  boundary,  y  =  1.  are  defined  as  reflecting  walls,  and  at  the  right 
boundary,  x  =  2.  outflow  boundary  conditions  are  imposed.  In  particular,  the  cells  at  the  exit 
plane,  at.  x  =  2.  are  assigned  values  for  the  flow  invariants  identical  to  those  just  upstream  of  the 
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end  cells.  The  interpolation  and  integration  schemes  using  the  Riemann  Invariant  Manifolds  can 
then  compute  the  fluid  fluxes  at  the  exit  plane  as  if  it  were  an  internal  interface.  This  procedure 
effectively  ignores  information  that  may  be  entering  the  domain  at  z  =  2  from  the  region  *  >  2. 
Inflow  boundary  conditions  are  specified  at  the  left  boundary,  x  =  0,  and  they  are  chosen  to 
match  the  experimental  conditions  as  closely  as  possible.  The  splitter  plate  separates  the  low- 
speed  (upper)  stream,  0.5  <  y  <  1,  designated  by  the  subscript  ‘1’,  from  the  (lower)  high-speed 
stream,  0  <  y  <  0.5,  designated  by  the  subscript  ‘2’. 

Two  such  runs  are  presented  in  the  following  figures.  Mach  number  contours  are  given  for  two 
cases.  In  both  cases  the  low  speed  side  corresponds  to  Mach  number  M\  =  0.65  and  the  high  speed 
side  to  Mach  number  M2  =  1.13.  The  high  speed  stream  enters  the  test  section  at  a  5  degree 
angle.  The  first  case,  shown  in  III.5,  corresponds  to  an  inlet  pressure  ratio  of  pi/p2  =  1.05  and  the 
second  case,  shown  in  III.6,  corresponds  to  a  pressure  ratio  of  p\/p2  =  1.5.  The  initial  condition  in 
the  entire  domain  is  that  of  the  inflow  conditions  and  the  computation  is  carried  out  with  a  CFL 
number  of  0.40  until  a  steady  state  is  achieved. 

1.0 
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.6 
.4 


.2 
.0 
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Fig.  III.5.  Mach  number  contour  levels  for  a  pressure  ratio  of  p\Jp2  =  1.05 .  There  are  60  contour 
levels  in  the  range  0.6  <  M  <  1.4.  The  resolution  is  100  X  50  cells. 


As  also  found  experimentally  (c/.  III.l,  III.  1),  the  main  feature  of  the  computed  flows  is  that 
the  subsonic  low-speed  stream  is  accelerated  and  becomes  supersonic.  The  shear  layer  curves  and 
acts  as  the  lower  boundary  of  a  gas-dynamic  nozzle  for  the  low-speed,  initially  subsonic,  side.  In 
the  first  case,  a  normal  shock  downstream  makes  the  flow  subsonic  again,  whereas  in  the  second 
case  an  oblique  shock  wave  system  is  set  up  and  the  low-speed  stream  remains  supersonic  through 
the  exit  of  the  computational  domain. 

III.4  Computation  of  chemically-reacting  flows 

As  part  of  our  effort  to  extend  these  computations  to  chemically-reacting,  compressible,  shear- 
layer  flows,  we  are  revisiting  the  original,  Lagrangian,  1-D  formulation  (Lappas  et  al.  1993),  explic¬ 
itly  incorporating  finite-rate  chemistry  in  the  system  equations.  At  this  writing,  we  have  developed 
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Fig.  III.6.  Mach  number  contour  levels  for  pressure  ratio  p\jp2  =  1-5.  60  contour  levels  in  the 
range  0.6  <  M  <  1.4.  The  resolution  is  100  x  50  cells. 

the  algorithm  for  the  computation  of  1-D  unsteady  flows  with  chemical  reactions  which  treats  the 
fluid  mechanics  and  the  chemistry  simulteneously  in  a  robust  way. 

Conventional  approaches  use  the  theory  of  characteristics  to  develop  a  number  of  numerical 
schemes  for  integrating  this  system.  In  such  schemes,  the  chemical  source  terms  are  then  integrated 
separately,  in  a  variety  of  empirical  ways.  Usually  they  are  treated  in  a  separate  step  and  con¬ 
sidered  as  a  “source”  term  for  the  non-reacting  flow.  The  theory  of  Riemann  Invariant  Manifolds, 
developed  in  the  previous  stages  of  this  project,  however,  shows  that  for  multidimensional  flow 
the  characteristic  surfaces  offer  no  advantage,  computational,  or  otherwise.  This  is  more  so  the 
case  in  reacting  flows.  Although  one  can  still  define  characteristic  surfaces,  they  no  longer  have  a 
simple  physical  interpretation,  because  one  cannot  define  a  speed  of  sound  in  flows  where  locally 
there  is  no  chemical  equilibrium.  Furthermore,  the  use  of  conventional  methods  can  not  handle 
the  additional  problem  of  the  very  large  range  of  time  and  length  scales  (computational  stiffness) 
associated  with  chemically-reacting  flows. 

The  new  algorithm  can  be  considered  as  an  extention  of  the  theory  of  Riemann  Invariant 
Manifolds.  The  key  idea  is  to  identify  trajectories  in  spacetime  on  which  appropriate  invariants 
exist  for  the  chemically-reacting  case.  This  allows  us  to  integrate  the  governing  equations  in  a 
mathematically  consistent  way.  We  expect  the  first  results  in  the  course  of  the  next  few  months. 

This  work  is  part  of  the  Ph.D.  research  effort  of  M.  Papalexandris. 
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CHAPTER  IV 

CHEMISTRY  IN  NONUNIFORM  FLOW 

IV.l  DISSOCIATION  RATES  WITH  VIBRATIONAL  NONEQUILIBRIUM 


Objectives  and  Status  of  Research 

One  of  the  primary  uncertainties  in  the  modeling  of  hypersonic  reacting  flows  is  how  the 
vibrational  state  of  a  molecule  affects  its  dissociation  rate.  Therefore,  the  main  objective  of  this 
research  is  to  test  chemical  reaction  rate  models  in  the  presence  vibrational  nonequilibrium.  To  do 
this,  we  are  using  experimental  results  from  the  Caltech  T5  shock  tunnel  to  assess  the  vibration- 
dissodation  coupling  models  that  are  used  in  computational  fluid  dynamics  codes.  These  models 
will  then  be  used  to  design  further  experiments  in  T5  to  refine  our  understanding  of  the  effect  of 
vibrational  nonequilibrium  on  reaction  rates.  These  new  experiments  will  determine  if  the  current 
models  are  lacking  in  their  ability  to  predict  important  flow  features.  If  not,  this  set  of  experiments 
will  provide  assistance  in  developing  new,  more  physically  consistent  models. 

To  date,  all  of  the  currently  used  vibration-dissodation  coupling  models  have  been  implemented 
in  a  two-dimensional  thermo-chemical  nonequilibrium  computational  fluid  dynamics  method.  The 
models  that  have  been  tested  are  the  Park  TTV  model,  the  CVD  and  CVDV  models  of  Treanor 
and  Marrone,  the  model  of  Macheret  and  Rich,  and  the  model  of  Olynick  and  Hassan.  The  flow  of 
reacting  nitrogen  over  a  sphere  has  been  simulated  using  these  models  and  compared  to  available 
experimental  data  from  T5.  The  shots  were  made  in  a  high  density  flow  regime  where  the  vibrational 
nonequilibrium  zone  is  small  compared  to  the  entire  flow  field.  Under  these  conditions,  it  was 
found  that  while  the  various  models  gave  different  vibrational  temperature  distributions,  the  shock 
stand-off  distances  and  density  contours  are  very  similar.  Thus,  these  comparisons  with  previous 
experiments  show  that  all  of  the  vibration-dissociation  models  are  able  to  predict  the  shock  stand¬ 
off  distances  and  density  fields  correctly.  Fig.  IV.l  shows  both  the  experimental  interferogram  and 
a  computed  interferogram  for  the  Park  model.  The  computed  shock  stand-off  distance  is  accurate 
to  within  a  tenth  of  a  millimeter,  while  the  fringe  patterns  are  virtually  identical. 

A  new  experiment  is  being  designed  for  T5  that  will  be  more  sensitive  to  the  differences  in  the 
reaction  rate  models.  Simulations  have  been  made  of  T5  flows  over  spheres  at  lower  densities,  and 
differences  among  the  models  are  discernible  in  both  the  shock  stand-off  distances  and  the  density 
contours.  Fig.  IV.2  shows  the  computed  density  contours  for  the  Park  and  the  CVDV  models  at  a 
density  of  10-5  kg/m3.  It  can  be  seen  that  the  shapes  of  the  contours  in  the  shoulder  region  differ 
dramatically.  Future  experimental  results  will  be  used  to  determine  which  of  the  models,  if  any, 
give  the  correct  results. 

In  addition,  computational  studies  are  being  carried  out  to  study  the  effect  of  the  choice  of 
oscillator  model  and  relaxation  time  expression  on  the  vibration-dissociation  coupling  problem. 
Once  appropriate  test  cases  have  been  determined,  experimental  results  will  be  able  to  determine 
which  choices  yield  the  best  results. 


Fig.  IV. 1.  Comparison  of  experimental  (left)  and  computational  (right)  differential  interferograms 
of  Nitrogen  flow  over  a  sphere.  Free-stream  conditions:  =  5070 m/s,  = 

0.0402  kg/m3,  Too  =  2260  K. 
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Fig.  IV.2.  Comparison  of  the  CVDV  model  (left)  and  the  Park  model  (right)  for  the  same  conditions 
as  those  of  IV. 1  except  /)<*,  =  0-0402  kg/m3. 
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IV. 2  ELECTRON-DRIVEN  REACTION  RATES  IN  HYPERSONIC  FLOW 


Background  and  Objective 

Electrons  resulting  from  ionization  due  to  bow-shock  heating  play  an  important  role  in  deter¬ 
mining  the  physical  and  chemical  properties  of  hypersonic  flows. 

Processes  of  interest  include  vibrational  excitation,  formation  of  negative  ions  (which  may  then 
fragment),  and  electron  impact  excitation  and  deexdtation  of  excited  metastable  states.  Further¬ 
more,  with  temperatures  as  high  as  10,000  Kelvin  due  to  non-equilibrium  overshoots,  electrons  in 
the  tail  of  the  energy  distribution  function  can  lead  to  electronic  excitation  of  ground  state  species, 
particularly  those  that  are  vibrationaliy  excited.  The  available  data  base  for  these  electron-molecule 
collision  processes  which  can  occur  in  flowfields  is  generally  fragmentary,  and,  in  many  cases,  in¬ 
consistent.  Furthermore,  cross  sections  for  electronic  excitation  of  molecules  by  electrons  with 
barely  enough  energy  to  drive  the  excitation  and  for  electron  collisions  with  electronically  excited 
molecules  and  with  radicals  such  as  OH  are  essentially  non-existent  and  experimentally  inaccessible. 
Finally,  the  experimental  effort  to  measure  electron-molecule  collision  cross  sections  is  very  limited 
and  shrinking  in  the  United  States. 

An  obvious  strategy  for  obtaining  these  electron  collision  cross  sections  is  to  calculate  them. 
However,  the  significant  computational  demands  of  such  calculations  have  made  such  an  approach 
impractical  to  date  with  conventional  (sequential)  computers.  The  objective  of  our  work  is  to 
exploit  the  high-performance  computing  provided  by  parallel  computers,  consisting  of  arrays  of  com¬ 
mercially  available  microprocessors,  to  calculate  the  cross  sections  for  electron-molecule  collisions 
needed  in  modelling  the  chemical  and  physical  properties  of  hypervelocity  flows.  With  appropriate 
computational  strategies,  the  high  aggregate  speeds  and  large  memory  of  these  parallel  computers 
make  such  an  approach  feasible. 


Progress  and  Highlights 

At  the  relatively  low  impact  energies  of  interest  in  hypersonic  flowfields,  the  accurate  treatment 
of  the  electron-molecule  collision  process  is  a  computationally  challenging  problem.  In  these  studies 
we  employ  a  multichannel  extension  of  the  variational  principle  for  scattering  amplitudes  originally 
introduced  by  Schwinger.1  This  formulation,  known  as  the  Schwinger  multichannel  method,2  was 
specifically  designed  to  address  low-energy  electron-molecule  collisions. 

As  in  most  variational  approaches,  there  are  two  major  computational  tasks  in  these  studies: 
(»)  calculation  of  the  primitive  matrix  elements  of  the  operators  arising  in  the  variational  expression 
for  a  given  expansion  basis,  and  (ii)  transformation  of  these  primitive  matrix  elements  to  gener¬ 
ate  the  physical  matrix  elements  of  the  linear  equations  defining  the  variational  solution.  In  the 
Schwinger  multichannel  method,  as  in  the  original  Schwinger  method,  the  trial  wave  function  need 
not  satisfy  scattering  boundary  conditions;  square-integrable  functions  such  as  Cartesian  Gaussians 
may  be  employed.  In  this  basis  all  matrix  elements  in  the  variational  expression  can  be  evaluated 
analytically,  except  those  involving  the  Green’s  function.  These  Green’s  function  matrix  elements 
have  no  known  analytic  form  and  must  be  done  by  quadrature.  In  fact,  evaluation  of  these  matrix 
elements  becomes  the  computationally  intensive  step  of  our  calculations.  The  occurrence  of  these 
difficult  matrix  elements  in  variational  principles  based  on  integral  equations  has  discouraged  their 
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use  to  date,  in  spite  of  highly  desirable  features  of  these  approaches  such  as  the  flexibility  of  em¬ 
ploying  trial  functions  which  do  not  satisfy  scattering  boundary  conditions  and  their  very  favorable 
convergence. 

Though  initially  implemented  on  conventional  (sequential)  computers,  the  computationally 
intensive  step  of  these  calculations  is  ideally  suited  to  computers  of  distributed-memory  parallel 
architectures.  Its  most  compute-intensive  tasks  are  (i)  the  evaluation  of  a  very  large  number 
(typically  10*°  or  more)  of  elementary  integrals  describing  the  interaction  of  pairs  of  electrons,  a 
step  requiring  a  code  of  about  2000  lines  of  Fortran,  and  (»»)  the  transformation  of  these  elementary 
quantities  into  the  matrix  elements  appearing  in  the  variational  expression.  The  evaluation  of  the 
two-electron  integrals  is  readily  parallelized:  one  simply  has  each  processor  independently  calculate 
a  different  subset  of  the  integrals.  Step  (tt),  the  transformation  of  the  integrals,  can  be  organized 
as  a  series  of  large  distributed  matrix  multiplications  between  the  global  array  of  integrals  and  a 
distributed  transformation  matrix  that  is  easy  to  construct. 

With  this  strategy  for  parallelization  of  our  computation,  we  have  achieved  sustained  per¬ 
formance  of  about  4  GFLOP  on  the  full  Intel  Delta,  a  parallel  computer  consisting  of  over  500 
i860  microprocessors.  This  is  approximately  100  times  the  performance  of  the  original  sequen¬ 
tial  program  on  the  CRAY  Y-MP.  Recently  we  have  successfully  ported  this  parallel  code  to  the 
CRAY  T3D,  a  multicomputer  with  128  DEC  alpha  processors.  On  this  machine  we  have  achieved 
performance  approximately  50  times  that  of  the  CRAY  Y-MP. 

We  have  exploited  this  computational  procedure  and  the  parallel  computer  resources  available 
at  the  Institute  and  the  Jet  Propulsion  Laboratory  to  carry  out  the  following  calculations  relevant 
to  the  modelling  of  hypersonic  flows: 

•  We  have  studied  the  cross  sections  for  near-threshold  excitation  of  the  long-lived  A  3  state 
of  N2  by  electrons.  These  near-threshold  excitation  processes  are  difficult  to  study  experimentally 
due  to  the  very  low  post-collision  energy  of  the  electron.  For  this  system  some  limited  cross  section 
data  based  on  detection  of  the  metastable  state  itself  is  available  at  about  2  eVs  above  threshold.3 
Comparison  of  our  calculated  cross  sections  with  this  data  is  very  encouraging. 

•  We  have  completed  the  first  calculation  of  the  cross  sections  for  excitation  of  the  B3nff  state 
of  N2  by  electron  impact  on  its  long-lived  A  3  state.  We  believe  that  such  cross  sections  for 
electron  impact  on  electronically  excited  metastable  states  can  play  an  important  role  in  modelling 
the  chemical  and  physical  properties  of  hypersonic  flows.  These  cross  sections  would  be  generally 
inaccessible  experimentally.  Calculations  of  these  cross  sections  required  significant  changes  in  our 
computer  codes  and  additional  calculations  of  these  cross  sections  are  under  way. 

•  We  have  obtained  preliminary  results  for  near-threshold  excitation  of  the  lowest  excited 
electronic  state  (3  £)*)  of  CO2. 

In  summary,  we  believe  that  these  applications  of  highly  parallel  computers  to  the  calculation 
of  electron-molecule  collision  cross  sections  are  an  early  indication  of  the  revolutionary  impact  that 
parallel  computers  can  be  expected  to  have  on  our  ability  to  robustly  model  complex  physical  and 
chemical  systems. 
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Publication 

The  results  of  our  calculations  of  the  electron  impact  excitation  cross  sections  for  No  and  CO2 
outlined  above  will  be  prepared  for  publication  within  the  next  six  months.  We  expect  these  results 
to  lead  to  two  publications. 

Personnel  associated  with  the  Research 
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volved  in  this  research  effort. 
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The  research  outlined  above  has  been  presented  at  the  following  seminars  and  meetings: 

•  Colloquium  speaker  at  the  Institute  for  Space  and  Astronautical  Sciences,  Tokyo,  Japan 
(July,  1993) 

•  Colloquium  speaker  at  the  Institute  for  Physical  and  Chemical  Research  (RIKEN),  Tokyo, 
Japan  (July,  1993) 

•  Invited  speaker  at  the  Gaseous  Electronics  Conference,  Montreal,  Canada  (October,  1993) 
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IV.S  NONEQUILIBRIUM  LEEWARD  SHOCK- VORTEX  AERODYNAMICS 


Objectives  and  Status  of  Research 

The  objectives  of  this  research  are  to  study  the  hypervelocity  flow  at  enthalpy  sufficient  to 
activate  dissociation/recombination  chemistry  about  cones  and  slender  bodies  at  large  incidence, 
with  particular  reference  to  the  effect  of  chemistry  on  the  structure  of  the  leeward  flow  shock/vortex 
system.  Significant  progress  has  been  made  in  two  relevant  areas.  These  are 

(1)  Completion  of  the  detailed  study  of  the  pseudo-Euler  flow  past  a  cone  at  incidence.  This  is 
described  in  detail  in  a  forth-coming  publication  (Macrossan  &  Pullin,  1994). 

(2)  The  development  and  testing  of  a  flexible,  viscous,  three-dimensional  (3D)  frozen  flow  code, 
PGPSD ,  on  the  Caltech  Intel  Paragon  parallel  computer.  Detailed  testing  of  this  code  is 
described  by  (Mallett,  Macrossan  &  Pullin,  1994)  currently  submitted  to  J.  Comp.  Phys. 
while  calculation  of  the  frozen  inviscid  and  viscous  flow  about  a  blunt  nosed  delta  wing  at 
incidence  is  given  in  (Mallett,  Pullin  &  Macrossan,  1994),  submitted  to  the  AIAA  J.. 

The  next  stage  of  the  work,  namely  the  incorporation  of  the  Lighthill-FYeeman  ideal-dissociating 
gas  chemistry  model  into  PGPSD  is  planned  for  the  near  future.  A  comprehensive  statement  of 
significant  accomplishments  in  the  areas  (1)  and  (2)  above  is  outlined  below.  These  statements 
constitute  the  current  status  of  the  research  effort.  They  are  followed  by  a  brief  statement  of  the 
expected  progress  within  the  next  year. 

Personnel  associated  with  the  Research 

(1)  D.  I.  Pullin,  Professor  of  Aeronautics. 

(2)  Dr.  M.N.  Macrossan,  University  of  Queensland,  Australia.  Collaborator. 

(3)  Dr.  R.  Mallett,  Postdoctral  Fellow  1993.  Consultant.  Author  of  principal  3-D  inviscid  and 
viscous  code. 

(4)  Mr  Shaun  Shariff,  Postgraduate  student.  Phd  project  will  probably  be  the  viscous  cone  flow 
with  reaction  chemistry. 

IV.3.1  Inviscid  Cone  Flow  with  IDG  Chemistry 

This  study  was  intended  as  a  precursor  of  the  more  complicated  viscous  flows  about  cones 
at  incidence  and  about  a  blunt-nosed  delta  wing.  The  computations  were  performed  with  the 
IDG  model  of  reacting  Nitrogen  about  a  15°  half-angle  cone  at  30°  incidence.  Flow  conditions 
were  Moo  =  5.27  Uoo  =  6.13  km/sec.  Several  runs  were  done  each  with  different  free-stream 
dissociation  fraction  a^,  density^,  in  the  range  poo  =  1  —  10-3fcpm-3  and  free-stream  enthalpy 
Ho  in  the  range  22  —  28  MJ  kg~i,  chosen  to  generate  a  range  of  dissociation  chemical  length  scales 
Ls  =  Uoo/(da/dt)9  where  “s”  refers  to  post- windward  shock  conditions.  Each  run  was  performed 
over  a  very  long  cone  in  order  to  study  the  transition  from  frozen  flow,  at  axial  lengths  x  «  Ls 
where  the  flow  will  be  nearly  frozen  (no  active  chemistry)  through  chemical  nonequilibrium  to 
equilibrium  flow  for  x  >>  Ls. 

Detailed  analysis  of  the  flow  in  the  windward  shock  layer  and  near  the  leeward  shock-vortex 
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system  have  been  performed  and  many  flow  diagnostics,  such  as  the  variation  following  a  particle 
path  along  a  streamline,  of  the  dissociation  fraction  a,  and  of  the  individual  dissociation  and 
recombination  reaction  rates,  have  been  studied,  and  are  reported  in  detail  in  Macrossan  &  Puliin 
(1994).  We  select  for  special  comment  here  several  flow  features  which  we  plan  to  investigate 
in  detail  for  the  viscous  case.  The  first  is  the  relative  stand-off  position  of  the  windward  shock 
from  the  cone.  It  was  found  that  the  quantity  Aa/z)/(p4/poo)>  where  (A,  is  the  standoff  distance 
scaled  against  distance  from  the  apex  z  and  pa  poo  are  the  free-stream  and  post-shock  densities 
respectively),  scales  well  against  z/I>4.  This  scaling  is  based  on  that  used  by  Hornung  ( J .  Fluid 
Mech.  53,  149,1972)  in  correlating  the  shock  standoff  distance  for  hypervelocity  reacting  flow  over 
spheres.  It  gives  a  good  correlation  up  to  xf  La  ss  2  at  which  point  the  data  for  different  stagnation 
enthalpies  separate,  leading  to  an  approach  to  the  equilibrium  limit  which  depends  on  Hq.  This  has 
been  analysed  by  Wen  (Phd  thesis,  California  Institute  of  Technology,  1994)  for  flow  over  spheres. 
A  related  analysis  for  the  cone  flow  should  be  possible  and  is  planned. 

Second,  we  discuss  the  dimensions  of  the  leeward  vortex,  the  lateral  scale  of  which  was  found  to 
vary  significantly  with  the  degree  of  active  chemistry.  When  plotted  as  Av/x 2  versus  x/L,  where 
Av  is  a  measure  of  the  vortex  area  in  a  plane  normal  to  the  cone  axis,  good  collapse  for  near  frozen 
flow,  and  up  to  x/La  «  50  was  found,  followed  by  a  separation  towards  the  equilibrium  limit. 

A  general  conclusion  for  the  inviscid  cone  flow  is  that  Ls  is  an  important  scaling  parameter  for 
the  effects  of  chemistry  in  both  the  windward  and  the  leeward  flow.  However  as  the  body  length 
becomes  large  compared  with  La ,  this  quantity  is  not  the  dominant  parameter.  A  one-dimesional 
streamtube  model  has  been  used  to  determine  a  more  appropriate  length  scale  Lr  j,  for  the  leeward 
flow,  when  the  equilibrium  limit  is  approached. 

IV.3.2  Viscous  3-D  Code  PGP3D 

A  code  for  the  computation  of  3-D  hypersonic  flows  ( PGP3D:  Parallel  General  Purpose  3D)  has 
been  implemented  on  the  Caltech  Intel  Paragon  computer.  The  code  has  been  written  for  maximum 
flexibility.  It  is  structured  as  a  generalized  finite  volume,  explicit  time-marching  method.  PGPSD 
currently  contains  switches  to  select  either  either  pseudo-inviscid  of  viscous  laminar  flow.  The 
flux-function  routines  are  structured  in  such  a  way  that  any  numerical  flux  function  for  computing 
fluxes  across  cell  interfaces  can  be  utilized  by  incorporating  appropriate  subroutines.  At  present 
two  choices  are  available,  either  the  equilibrium  flux  method,  a  kinetic  based  scheme  based  on 
flux-split  at  the  particle  distribution  function  level,  or  a  fast  Riemann  solver  based  in  part  on  the 
Roe  algorithm,  with  provision  for  a  full  Godunov  treatment  of  strong  shocks.  These  methods  have 
been  tested  in  detail  against  standard  results  for  the  hypersonic  flat-plate  boundary  layer  problem. 
This  work  is  described  in  a  publication  (Mallett,  Macrossan  &  Puliin,  1994)  currently  submitted 
to  J.  Comp.  Phys. 

Detailed  computations  have  been  done  of  the  pseudo-inviscid  and  viscous  flow  about  a  blunt 
nosed  delta  wing  of  70°  sweep,  at  an  angle  of  attack  of  30°.  The  freestream  Mach  number  was 
fixed  at  8.7.  For  the  laminar  viscous  flow  two  Reynolds  numbers  of  Re  —  2.25  X  105  and  2.25  x  I06, 
with  the  initial  aim  of  comparing  pseudo-inviscid  and  viscous  leeward  flow  at  large  Re.  This  work 
is  described  in  detail  in  a  publication  (Mallett,  Puliin  &  Macrossan,  1994)  currently  under  review 
by  the  A1AA  J..  Fig.  refxl  shows  a  comparison  of  Mach  number  contours  for  the  pseudo-inviscid 
flow  and  the  two  viscous  flows  at  a  station  along  the  delta-wing.  For  the  inviscid  flow,  a  strong  and 
curved  cross-flow  shock  is  apparent  on  the  leeward  side,  whilst  for  the  viscous  flow,  the  leeward 
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Fig.  IV.3.  Contours  of  Mach  number  at  z/L  =  0.5.  Left;  pseudo-inviscid.  Middle;  Re  =  2.25  x  10®. 
Rjght;Re  =  2.25  x  10® 

shock  is  merged  with  a  separating  shear  layer. 

Future  work 

In  order  to  further  test  PGPSD,  it  is  planned  in  the  near  future  to  repeat  at  least  the  Re  = 
2.25  x  10®  laminar  viscous  flow  calculation  using  the  Fast  Riemann  solver  in  place  of  the  EFM 
numerical  flux  function,  and  to  make  detailed  comparisons  of  the  same  flow  as  calculated  by  the 
two  methods.  At  the  same  time  the  IDG  chemistry  model  will  be  incorporated  into  PGPSD  in  order 
to  investigate  the  effect  of  chemistry  on  the  windward  and  leeward  flows  for  both  the  cone-flow  and 
the  delta-wing  flow  with  particular  reference  to  the  shock  stand-off  distance,  the  heat  transfer  on 
the  windward  body  surface,  and  the  structure  of  the  leeward  flow.  In  the  longer  term  it  is  planned 
to  incorporate  a  more  complex  partial  equilibrium  model  of  air  (Rein,  Phys.  Fluids  ,  A,  4,  873 
1992)  into  PGPSD.  Shaun  Shariff  is  currently  gaining  experience  with  a  2D  version  of  the  code 
PGPSD  and  it  is  planned  that  he  will  move  on  to  the  viscous  reacting  cone  flow  in  1995. 

Publications  resulting  from  project 

Macrossan,  M.N  and  Pullin,  D.I.,  A  computational  investigation  of  inviscid,  hypervelocity  flow  of 
a  dissociating  gas  past  a  cone  at  incidence.  Accepted  for  publication  in  J.F.M.  In  Press,  1994. 

Mallett,  E.R.  A  numerical  study  of  hypersonic  flow  over  a  blunt-nosed  delta  wing  at  incidence.  Phd 
Thesis,  University  of  Queensland,  Jan  1994. 

Mallett,  E.R.  Macrossan,  M.N.  and  Pullin,  D.I.,  The  equilibrium-flux  method  for  the  solution  of 
the  Navier-Stokes  equations  for  a  compressible  ideal  gas.  Submitted  to  J.  Comp.  Phys. 

Mallett,  E.R.  Pullin,  D.I.  and  Macrossan,  M.N.,  A  numerical  study  of  hypersonic  leeward  flow  over 
a  blunt  nosed  delta  wing.  Submitted  to  AIAA  J. 
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CHAPTER  V 

DIAGNOSTICS 


V.l  Introduction 

Since  the  major  part  of  the  research  in  this  program  is  experimental,  much  of  our  work  is  also 
oriented  toward  development  of  better  diagnostic  methods.  In  some  instances  the  techniques  bang 
developed  are  very  closely  coupled  to  particular  projects  and  the  progress  in  their  development 
can  be  more  effectively  presented  together  with  the  progress  in  the  sub-project  in  which  they  are 
being  applied.  This  is  the  case  for  the  development  of  the  new  thermocouple  heat  transfer  gauges 
described  in  section  II.2,  and  the  new  high-signal-to-noise  laser  scattering  technique  described  in 
section  III. 

In  the  case  of  the  following  section  on  laser-induced  thermal  acoustics,  the  presentation  is  separated 
from  a  specific  application,  because  it  needed  to  be  developed  as  a  project  in  its  own  right  at  firs. 
The  method  is  now  maturing,  however,  and  will  be  applied  to  the  other  sub-projects  near  the  end 
of  the  next  year  of  this  program. 

V.2  DIAGNOSTICS  WITH  LASER-INDUCED  THERMAL  ACOUSTICS 

Objectives  and  Status  of  Research 

The  objectives  of  the  research  described  in  this  subsection  are  to  develop  laser-induced  thermal 
acoustics  (LITA),  a  nonlinear  optical  technique  for  single-shot  measurement  of  gas  and  flow  prop¬ 
erties,  including  species  concentration,  temperature,  transport  properties,  and  velocity,  and  to  apply 
LITA  to  hypersonic  flows  in  T5.  These  objectives  are  extended  from  those  in  the  URI  proposal  as 
a  consequence  of  favorable  developments  during  LITA  research.  The  feasibility  of  this  technique 
and  the  accuracy  of  an  expression  derived  for  the  LITA  signal  (Ref.  1)  were  verified  in  a  series  of 
experiments  on  laboratory  air.  A  second  series  of  experiments,  currently  under  construction,  will 
develop  LITA  concentration  and  velocity  measurements  using  a  small  shock  tunnel  and  blow-down 
supersonic  wind  tunnel.  These  experiments  are  also  for  preparation  and  risk-reduction  for  LITA 
studies  in  T5,  which  are  planned  for  early  1995. 

Equipment  for  LITA  studies  furnished  through  URI  and  Caltech:  narrowband  dye  laser,  Cynosure 
flashlamp-pumped  dye  laser  and  modification,  high-pressure  test  section  with  optical  access,  high- 
resolution  (<  3  GHz)  echelle  spectrometer,  Princeton  Instruments  thermo-  electrically  cooled  CCD 
camera  and  recording  system,  Tektronix  four-channel  2  GSamples/s  digital  storage  oscilloscope, 
Wideband  amplifed  PMT  optical  detection  system,  optics  and  mounts,  triggering  electronics. 

V.2.1  Description  of  the  LITA  technique 

Over  the  past  year,  work  on  optical  diagnostics  for  T5  has  progressed  markedly  with  the  develop¬ 
ment  at  GALCIT  of  LITA,  a  new  technique  superior  in  many  respects  to  “conventional”  techniques. 
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In  the  first  series  of  experiments  LITA  has  provided  relatively  simple,  accurate,  single-shot,  spa¬ 
tially  resolved  measurements  of  sound  speed,  transport  properties,  and  collisional  thermalization 
rate  in  laboratory  air.  Experiments  curr  :,!'v  under  development  will  apply  LITA  for  single-shot 
measurement  of  these  properties,  as  well  as  ‘;.e  complex  gas  susceptibility,  susceptibility  spectra, 
and  the  gas  velocity.  Absolute  species  concentration  and  temperature  measurements  may  be  ex¬ 
tracted  from  susceptibility  spectra.  The  work  on  LITA  performed  at  GALCIT  has  the  interest  of 
the  wider  combustion  di agnostic  community  because 

•  quantitative  signal  analysis  is  simple  (cf.,  LIF,  DFWM), 

•  exhaustive  information  about  probed  species  is  not  required  (cf.,  LIF), 

•  LITA  excels  in  strongly  quenching  environments  (e.g.,  where  collision  rates  are  high,  a 
requirement  for  fast  chemistry;  cf.,  LIF), 

•  LITA  requires  only  standard  lasers  and  equipment,  and 

•  LITA’s  capabilities  exceed  those  of  other  techniques  (e.g.,  single-shot  transport-property, 
absolute  concentration,  and  velocity  measurements) 

•  LITA  signals  are  relatively  strong. 

LITA  is  a  four- wave  mixing  (FWM)  technique  related  to  coherent  anti-Stokes  Raman  scattering 
(CARS)  and  degenerate  FWM  (DFWM).  In  FWM  techniques,  crossed-laser  beams  form  electric- 
field  interference  “gratings”  in  the  probed  volume.  Electric  field  gratings  generate  susceptibility 
(index-of-refraction)  gratings  via  nonlinearities  of  the  probed  medium.  These  gratings  scatter  light 
from  a  third  beam  into  the  signal  beam.  If  the  intensity  of  the  signal  beam  can  be  related  to 
properties  of  the  medium,  FWM  can  be  used  to  measure  these  properties  remotely.  The  princi¬ 
pal  difference  between  the  techniques  is  the  type  of  nonlinearity  used  to  form  the  susceptibility 
gratings.  CARS  and  DFWM  form  gratings  by  quantum  mechanical  effects.*  These  gratings  have 
short  lifetimes  (<  5  ns  at  STP),  since  they  are  destroyed  by  molecular  collisions.  LITA  gratings  are 
acoustic:  crossed  driver  lasers  generate  acoustic  pressure  and  temperature  gratings  via  electrostric- 
tion  and  thermalization,  two  semi-classical  opto-acoustic  effects.  Acoustic  gratings  are  damped  by 
thermal  and  viscous  diffusion  on  a  much  longer  time  scale  (e.g.  1  /us)  than  quantum  mechanical 
gratings.  Longer  time  scales  permit  single-shot  time  resolution  of  the  evolution  of  the  signal  which 
allows  sound  speed  measurement.  Furthermore,  the  formation  and  lifetime  of  acoustic  gratings  are 
enhanced  by  quenching,  which  has  the  opposite  effect  on  quantum-mechanical  gratings. 

The  mathematics  of  acoustic  grating  formation  is  generally  far  simpler  than  that  of  quantum- 
mechanical  gratings.  This  makes  the  expression  relating  experimental  LITA  signals  to  gas  properties 
simpler  than  that  of  DFWM  or  CARS.  Thus,  it  was  possible  to  predict  the  LITA  signal  with 
an  analytical  expression  derived  previously  from  the  linearized  equations  of  light  scattering  and 
hydrodynamics.  By  fitting  the  experimental  signal  to  the  prediction,  sound  speeds  could  be  obtained 
to  an  accuracy  within  the  uncertainty  of  independent  determinations  (~0.5%).  Transport  properties 
obtained  this  way  matched  published  quantities  to  within  ~30%.  A  finite- beam-size  effect  neglected 
in  the  predictive  analysis  limited  the  accuracies  of  these  measurements.  The  analysis  has  since  been 


*  Thermal  or  acoustic  gratings  often  contribute  to  DFWM  signals,  especially  at  high  densities.  However,  the 
contribution  to  the  signal  from  these  gratings  is  generally  regarded  as  “interference”  with  the  (desired)  signals 
from  the  quantum-mechanical  gratings.  For  this  reason,  the  distinction  is  made  that  DFWM  gratings  are  purely 
quantum-mechanical  and  LITA  gratings  are  purely  acoustic,  i.e.,  LITA  signals  often  interfere  with  DFWM  signals. 
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extended  to  include  this  additional  effect,  the  effect  of  electrostriction,  and  the  effect  of  finite  laser 
linewidth.  Work  is  under  way  to  publish  the  extended  analysis  and  automate  LITA  data  reduction. 
This  maturity  in  signal  analysis  was  not  expected  until  late  in  the  URI  program. 


Fig.  V.l.  Schematic  diagram  of  LITA  operation.  For  notation  see  text. 

V.2.2  Results  of  demonstration  experiments 

Experiments  conducted  to  date  at  GALCIT  use  the  setup  schematized  in  Fig.  V.l.  A  pulsed 
narrowband  dye  laser  pumped  by  a  Q-s witched,  frequency-doubled  Nd:YAG  laser  drives  the  acoustic 
waves.  The  dye  laser  emits  about  40  mJ  of  light  tuned  between  587  nm  and  592  nm  in  approximately 
5  ns  pulses.  Three-millimeter-wide  driver  beams  are  formed  from  the  dye  laser  emission  by  a  50% 
beamsplitter  ( B.S .)  and  a  mirror  mj.  The  driver  beams  intersect  at  the  focus  of  the  400  mm  lens 
where  they  approximate  finite  plane  waves  about  200  pm  in  diameter.  The  distance  between  the 
parallel  beams,  adjusted  using  micrometer  translation  stages,  sets  the  beam-crossing  angle.  Three- 
dimensional  phasematching  (BOXCARS)  geometries  with  crossing  angles  ranging  from  1.0°  to  2.3° 
were  used  in  these  experiments.  The  acoustic  gratings  are  formed  in  the  sample  volume  defined  by 
the  intersection  of  these  beams.  Detection  of  the  driver  beam  pulse  by  a  silicon  PIN  photodiode, 
Pli  triggers  the  acquisition  of  the  LITA  signal.  Light  from  a  1W  CW  argon  ion  “source”  laser 
operating  at  488  nm  scatters  coherently  off  the  acoustic  gratings  into  the  signal  beam.  A  two-lens 
system  adjusts  the  source  beam  diameter,  which  ranged  from  about  200  pm  to  2  mm.  Phase¬ 
matching  adjustments  are  made  using  a  micrometer-driven  translation  stage  under  m*.  A  silicon 
PIN  photodiode,  P2,  monitors  the  intensity  of  the  source  beam. 

The  signal  beam  detection  system  consists  of  a  fast  photomultiplier  tube  (PMT)  with  at  least  a 
500  MHz  signal  bandwidth.  Optical  filters,  including  a  10  nm  interference  filter  at  488  nm  and  a  40 
mm  lens/20  pm  pinhole  spatial  filter,  prevent  signal  contamination.  The  LITA  signal  path-length 
to  the  detector  is  about  2.5  m.  Data  are  recorded  using  a  fast  digital  storage  oscilloscope  (DSO). 
This  records  typically  500-to- 2000-sample  time  histories  of  signals  from  the  PMT  and  P2- 

The  first  LITA  experiments  used  laboratory  air.  Nonresonant  electrostriction  generated  LITA 
signals  with  a  signal- to-source  beam  intensity  ratio,  or  reflectivity,  of  ~  10-8.  An  average  over 
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Fig.  V.2.  LITA  signals  averaged  over  256  shots  in  room  air  near  a  weak  line  of  NO2.  which  is  present 
to  about  50  ppb.  Full  line:  Signal  recorded  near  the  peak  of  the  spectral  line.  Dashed 
line:  Resonantly  enhanced  signal  near  the  spectral  line.  Dotted  line:  Non-resonant  signal. 
By  fitting  these  to  the  theoretical  prediction,  the  speed  of  sound  is  obtained  to  ±0.5%, 
and  the  bulk  viscosity  to  ±10%. 

256  shots  yields  the  nonresonant  data  plotted  as  the  dotted  curve  in  Fig.  V.2  .  The  modulation  or 
“ringing”  results  from  beating  between  Doppler-shifted  light  scattered  off  the  laser-induced  acoustic 
waves.  This  occurs  at  a  frequency  proportional  to  the  sound  speed  and  known  geometrical  and 
optical  parameters.  The  frequency  of  the  modulation  (about  26  MHz)  yields  a  sound  speed  accurate 
within  the  uncertainty  of  the  knowledge  of  the  laboratory  sound  speed  (~0.5%).  Electrostriction  is 
resonantly  enhanced  near  absorption  lines.  The  resonantly-enhanced  signal  recorded  near  a  weak 
line  of  NO2  appears  as  a  dashed  curve  in  Fig.  V.2  for  comparison.  When  the  driver  laser  is  tuned 
to  an  absorption  line  of  the  gas,  molecules  absorb  light  energy  and  release  thermal  energy  through 
inelastic  collisions,  generating  a  thermal  grating.  If  thermal  gratings  form  rapidly,  expansion  of 
the  gas  emits  acoustic  waves.  The  signal  recorded  near  the  peak  of  the  absorption  line  is  the 
solid  curve  plotted  in  Fig.  V.2,  which  shows  the  difference  in  signature  between  electrostriction 
and  thermalization  signals.  These  signals  are  proportional  to  the  real  and  imaginary  parts  of  the 
complex  gas  susceptibility  respectively,  allowing  measurement  of  the  complex  susceptibility.  LITA’s 
good  sensitivity  is  demonstrated  since  the  ambient  air  contained  less  than  50  ppb  of  NO2  according 
to  the  Southern  California  Air  Quality  Management  District. 

Seeding  the  gas  with  NO2  provided  LITA  signals  that  were  several  orders  of  magnitude  stronger 
than  the  signals  obtained  from  room  air,  with  estimated  reflectivities  around  10~4.  A  sample  single¬ 
shot  LITA  signal  appears  in  Fig.  V.3  with  the  theoretical  fit  plotted  with  a  dashed  curve,  showing 
excellent  agreement  with  experiment.  The  sound  speed  and  transport  coefficients  used  in  the  fit  are 
published  values  for  dry  air.  Here,  the  thermalization  process  is  treated  as  the  superposition  of  two 
single-rate  thermalizations:  a  fast  ~10  ns  process  and  a  slow  ~100  ns  process.  A  better  fit  to  the 
data  could  be  obtained  by  adopting  a  more  accurate  thermalization  model  and  source-beam-profile 
model  in  the  analysis  (a  top-hat  profile  is  assumed)  and  including  the  effect  of  source-beam  phase 
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Fig.  V.3.  Single-shot  LITA  record  in  room  air  seeded  with  NO2.  Dashed  curve  shows  theoretical 
prediction. 

noise.  These  improvements  are  being  incorporated  into  a  computer  fitting  routine  for  LITA  signals. 

Information  gained  from  these  experiments  has  guided  plans  for  future  experiments.  The  accuracy 
of  the  simple  analytical  model  and  the  relative  strength  of  LITA  signals  have  opened  new  paths  of 
exploration:  multiplex  LITA  and  LITA  velocimetry.  Absolute  species  population  and  concentration 
measurements  may  be  extracted  from  spectra  obtained  by  LITA  experiments  in  which  the  source 
laser  is  scanned  across  a  band  of  species  spectral  lines.  Such  scanning  techniques  are  common  in 
laser  diagnostics;  however,  they  are  not  amenable  to  pulsed-facilities  such  as  T5.  Alternatively,  a 
broadband  source  laser  can  be  used  along  with  dispersion  of  the  signal  beam  in  a  spectrometer  to 
obtain  single-shot  “multiplex”  spectra.  Multiplexing  using  conventional  FWM  techniques  has  two 
drawbacks  avoided  by  LITA. 

•  Analyzing  multiplex  spectra  from  “quantum-mechanical”  grating  techniques  is  complicated 
by  quantum  interferences,  etc.  Because  LITA  uses  classical  gratings,  quantum  mechanical 
complications  are  avoided.  The  simple  analysis  already  experimentally  verified  extends 
readily  to  multiplex  LITA. 

•  Much  of  the  energy  of  the  source  laser  falls  between  lines  and  is  wasted,  thus  a  high  reflec¬ 
tivity  is  needed.  Spectrometer  losses  compound  this  effect.  LITA  signals  with  reflectivities 
of  10-4  can  be  used  in  multiplex  experiments  studying  even  trace  species.  Furthermore, 
a  long-pulse,  broadband  source  laser  such  as  a  flashlamp-pumped  dye  laser  may  be  used. 
These  are  available  inexpensively  with  high  pulse  energies  (>  1  J). 

Multiplex  experiments  are  currently  under  construction.  A  flashlamp-pumped  dye  laser  with  about 
2.5  J  per  pulse  is  being  modified  for  less  mode  structure.  Laser  mode  structure  interferes  with  all 
multiplex  measurements.  A  high  resolution  echelle  spectrometer  with  multiple  inputs  has  been 
aligned  and  focused.  Spectra  fall  on  a  computer-controlled  precision  CCD  array.  The  spectrometer 
has  a  dynamic  range  of  at  least  2000:1  and  a  frequency  resolution  measured  to  be  better  than  3 
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GHz  at  589  nm  (>  105:1). 

LITA  velocimetry,  by  analogy  with  laser  Doppler  velocimetry  (LDV),  is  also  actively  being  pursued 
at  GALCIT.  In  this  technique,  the  Doppler  shift  of  light  scattered  off  convecting  acoustic  gratings 
is  measured  relative  to  an  unshifted  reference.  This  should  allow  single-shot  velocimetry  in  pulsed 
experiments  with  flow  velocities  from  50  to  3000+  m/s  (estimation  based  on  detector  frequency 
response,  beam  size  and  duration  effects  etc.)  For  these  experiments,  a  windowed  test  section  with 
optical  access  from  four  sides  has  been  designed  and  built.  The  test  section  is  designed  to  be  used  in 
connection  with  a  variety  of  test  articles  including  a  small  combustion-driven  shock  tunnel  capable 
of  generating  exit  velocities  of  3000  m/s  in  air  and  7000  m/s  in  hydrogen.  Alignment  and  testing 
will  first  be  performed  on  a  steady  ~Mach  2  jet  of  air  from  the  laboratory  air  supply. 

LITA  tests  in  T5  have  been  postponed  until  Spring  1995  to  allow  experimental  development  and 
exploration  of  the  potential  for  LITA  diagnostics.  The  risk  reduction  from  this  measured  devel¬ 
opment  and  knowledge  gained  about  the  use  of  this  promising  diagnostic  will  enhance  T5  LITA 
tests.  In  addition,  timely  development  of  generic  LITA  capabilities  such  as  spectral  multiplexing 
and  velocimetry  is  in  the  broad  interest  of  the  combustion  diagnostic  community. 
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